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Introduction 


^There  are  two  categories  of  magnetic  anomalies  associated  with  the  ocean.  Primary 
anomalies  are  induced  directly  by  interaction  of  ocean  flow  with  the  geomagnetic  field,  or 
induced  by  ionospheric  fluctuations,  or  are  the  expression  of  the  remnant  magnetization  of  the 
ocean  floor  rocks.  Secondary  anomalies  are  created  by  interaction  of  electric  current  flow  in  the 
ocean,  with  irregularities  of  tie  electrical  conductivity  distribution  within  the  ocean,  or  with 
irregularities  of  the  boundaries  of  the  ocean  (seafloor  ridges  and  troughs  or  ocean  surface 
waves).  Our  focus  in  this  report  is  on  secondary  anomalies,  although  some  attention  is  paid  to 
primary  effects  for  completeness.. 

Large  scale  and  relatively  intense  electric  currents  in  the  ocean  are  mainly  induced  by 
ionospheric  fluctuations.  The  primary  magnetic  fields  and  induced  electric  fields  in  the  ocean 
both  have  "reef  continuum  spectra  in  the  frequency  range  from  one  cycle  per  day  upward, 
which  is  caused  by  the  normal  unrest  of  the  ionosphere  and  by  bay  disturbances  and  the  princi¬ 
pal  part  of  magnetic  storms.  In  addition,  there  is  a  line  spectrum  consisting  of  harmonics  of  1 
cpd:  the  solar  daily  variation.  Typical  amplitudes  of  electric  fields  induced  in  the  ocean  are  of 
order  1  fiV/m,  although  larger  values  can  be  found  in  a  magnetic  storm. 

The  spatial  scale  of  the  anomalous  conductivity  structures  within  the  ocean  which  we 
shall  treat  are  relatively  small;  for  ocean  surface  waves  the  wavelength  is  typically  a  few.  hun¬ 
dred  meters,  and  a  comparable  scale  is  appropriate  for  ocean  structures  such  as  internal  waves 
and  frontal  systems.  Ocean  bottom  hills,  ridges,  and  troughs  are  typically  of  a  few  to  a  few  tens 
of  kilometers  in  size.  As  a  consequence  of  the  slow  variation  of  the  primary  electric  field  in  the 
jeean  and  the  small  spatial  scale  of  the  conductivity  anomalies,  the  magnetic  Reynold’s 
number,  R^ficrLD,  applicable  to  the  anomalous  secondary  fields,  is  small.  For  example,  if 
at,  the  oscillation  frequency,  is  10“4 see"1,  the  conductivity  of  seawater  or  is  3  5/m,  and  the 
horizontal  and  vertical  scales  L  and  D  are  40  km  and  4  km  respectively,  then  R-0.06.  From 
this  rough  calculation  it  is  apparent  that  time  variations  of  the  primary  fields  do  not  lead  to 
appreciable  self-inductance  effects  in  the  secondary  fields.  In  the  subsequent  calculations  we 
have  neglected  the  self-inductance  of  secondary  fields  entirely,  thus  treating  them  as  though 
they  were  caused  by  a  perfectly  static  primary  field.  This  approximation,  valid  for  these  secon¬ 
dary  fields,  is,  of  course,  entirely  inappropriate  for  the  treatment  of  the  primary  fields,  because 
their  large  horizontal  scale  ordinarily  leads  to  magnetic  Reynold’s  numbers  of  order  unity  or 
greater. 


Chapter  I 


Magnetic  Anomalies  Induced  by  Surface  Waves 


1.  Introduction 

Here  we  consider  magnetic  anomalies  associated  with  waves  of  small  amplitude  mov¬ 
ing  on  the  surface  of  a  deep  ocean.  There  are  two  principal  effects  in  this  category;  that  due  to 
interaction  of  the  surface  waves  with  a  uniform  electrical  current  flow  within  the  ocean  (such  as 
might  be  induced  by  electrical  activity  in  the  ionosphere),  and  that  due  to  wave-associated  sea¬ 
water  movement  through  the  Earth's  magnetic  held.  Section  2  of  this  chapter  deals  with  the 
magnetic  field  generated  by  the  uniform  current  flow  alone,  and  section  3  derives  the  perturb¬ 
ing  magnetic  effects  of  the  interaction  of  a  surface  wave  with  this  flow.  Section  4  is  a  note  in 
passing  that  the  mathematics  in  sections  2  and  3  is  also  applicable  to  the  calculation  of  magnetic 
anomalies  produced  by  a  special  class  of  bottom  features.  Section  5  discusses  the  interaction  of 
the  water  motions  that  accompany  surface  waves  with  the  Earth's  magnetic  field,  summarizing 
previous  work  by  others.  Section  6  covers  the  subject  of  the  gradients  of  the  magnetic  field 
vector  components  generated  by  the  mechanisms  explored  in  sections  2,  3,  and  S.  Finally,  sec¬ 
tion  7  puts  some  numbers  into  the  field  expressions  derived  in  sections  2,  3,  and  5  in  order  to 
give  some  idea  of  the  relative  magnitudes  of  the  different  magnetic  effects  under  various  cir¬ 
cumstances.  The  treatment  is  as  general  as  is  consistent  with  the  small  amplitude  wave,  deep 
uniform  ocean  model,  and  should  lend  itself  to  expansion  to  consider  the  effects  of  ocean  bot¬ 
tom  features  of  small  aspect,  should  this  be  desired. 
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2.  Uniform  horizontal  electric  field  In  a  flat-surfaced  ocean 

One  of  the  consequences  of  magnetic  field  fluctuations  within  the  volume  of  the 
ocean,  such  as  may  result  from  the  solar  daily  variation  of  the  Earth's  magnetic  field  or  from 
disturbances  in  the  upper  atmosphere,  is  the  generation  of  an  electric  potential  gradient  within 
the  ocean  and  the  corresponding  induction  of  an  electrical  current  flow;  for  example,  the  solar 
daily  variation  typically  generates  a  slowly  oscillating  electric  field  throughout  the  deep  sea  with 
an  amplitude  that  is  typically  a  large  fraction  of  a  microvolt  per  meter,  and  magnetic  storms  in 
the  upper  atmosphere  can  give  rise  to  electric  fields  several  times  larger.  Magnetic  fluctuations 
of  this  sort  are  for  the  most  part  of  a  sufficiently  ‘jw  frequency  that  the  magnetic  induction 
number  upal2  (where  a>  is  the  frequency  of  fluctuation,  n  and  or  are  the  magnetic  permeability 
and  electrical  conductivity  respectively  of  the  material  under  consideration,  in  this  case  seawa¬ 
ter,  and  1  is  the  length  scale  of  the  induction  effect  in  question,  in  this  case  the  depth  of  the 
ocean)  is  a  great  deal  less  than  unity,  indicating  that  the  magnetic  fluctuations  penetrate 
through  the  ocean  virtually  without  resistance;  as  a  result  of  this,  in  a  flat-bottomed  ocean  the 
induced  electric  field  will  be  essentially  uniform  with  depth.  Also,  and,  in  this  context,  perhaps 
more  to  the  point,  the  time  scale  of  these  magnetic  field  fluctuations  is  typically  sufficiently 
large  (with  periods  on  the  order  of  an  hour  and  up)  compared  to  that  of  water  wave  effects  on 
the  ocean  surface  that  the  electric  fields  induced  by  them  can,  for  our  present  purposes,  be  con¬ 
sidered  static  to  excellent  approximation.  Accordingly,  we  assume  a  horizontal,  spatially  uni¬ 
form,  temporally  static  electric  field  E<#  within  the  volume  of  the  ocean,  and  an  induced  electr¬ 
ical  current  density  distribution  J  such  that 

7-<r(z)£o*  ,  (2.1) 

where  <r(z)  is  the  electrical  conductivity  distribution  of  the  ocean,  assumed  here  to  be  a  func¬ 
tion  of  depth  only.  In  the  following  work  dealing  with  the  magnetic  effects  resulting  from  this 
induced  electrical  current  distribution,  deviation  of  the  current  density  from  /  as  given  by  (2.1) 
will  be  considered  as  due  to  the  perturbing  effects  of  variations  of  the  upper  and/or  lower  boun¬ 
daries  of  the  ocean,  and  for  a  flat-surfaced,  flat-bottomed  ocean  should  be  completely  absent. 

We  consider  first  the  magnetic  field  within  and  above  a  flat-surfaced,  flat-bottomed 
ocean  of  depth  D  resulting  from  the  induced  electrical  current  density  distribution  /  os  given 
by  (2.1).  If  this  current  distribution  is  viewed  as  the  limit  of  a  large  number  of  parallel  fine 
threads  of  electrical  current,  then  one  sees  that  the  resultant  magnetic  field  can  have  no  com¬ 
ponent  in  the  x  direction,  and  consideration  of  the  symmetry  of  the  problem  in  the  $>  direction 
shows  that  the  field  can  have  no  z  component,  where  z  points  vertically  upward,  This  field  can 
be  deduced  by  use  of  the  integral  form  of  Ampere's  law, 
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,  (2.2) 

together  with  some  symmetry  arguments,  and  the  total  magnetic  field  may  then  be  inferred  by 
adding  this  component  to  the  inducing  field,  which  can  be  derived  by  other  means. 

We  start  with  a  rectangular  path  of  integration  within  a  plane  perpendicular  to  £,  its 
top  and  bottom  legs  parallel  to  p  and  above  the  surface  and  below  the  bottom  respectively. 
From  the  above  arguments  we  have  that  the  induced  magnetic  induction  vector  B/  is  every¬ 
where  parallel  to  p,  and  from  considerations  of  symmetry  we  have  that  the  fields  above  the 
ocean  surface  and  beneath  the  bottom  should  have  the  same  magnitude  but  opposite  directions. 
Since  the  vertical  legs  of  the  rectangle  contribute  nothing  to  the  integral,  we  therefore  have  for 
the  magnetic  induction  everywhere  beneath  the  ocean  floor 

ti-Z—p  }<r(z)dz  ,  (2.3) 

*  o 

with  the  induced  magnetic  induction  everywhere  above  the  ocean  surface  being  just  the  nega¬ 
tive  of  this.  We  now  place  the  lower  leg  of  the  rectangle  exactly  at  the  ocean  floor,  and  move 
the  upper  leg  down  below  the  surface  to  a  depth  a  distance  z0  above  the  bottom;  given  (2.3)  for 
the  magnetic  induction  along  the  lower  leg,  we  have  from  (2.2)  the  magnetic  induction  along 
the  upper  leg 


s^—af-y 


f  tr(z)  dz  —  ?  w(z)  dz 
o  *0 


(2.4) 


If  the  electrical  conductivity  of  the  seawater  is  constant  over  the  depth  of  the  ocean  with  a 
value  <to,  then  this  reduces  to 

5i(z)  ~  —l/i  hqJqDP  ,  z>D  (2.5) 

5?  (2)  -  +V6  mo  h  C£»~2z0)  p  ,  0  <  z  <  D  (2.6) 

S',(z)-+'/2MJ()DP  ,  z  <  0  (2.7) 

where 


(neglecting  displacement  currents,  a  valid  approximation  in  this  context)  confirms  this  work, 

For  the  total  magnetic  induction  Br  at  the  ocean  floor  we  have 

BT-^y  ,  (2.10) 

where  Ze  is  the  E  mode  impedance  of  the  Earth,  a  convex  function  of  frequency  that  is 
dependent  in  ideal  circumstances  only  on  the  distribution  of  electrical  conductivity  below  the 
ocean  floor  (see  for  example  Cox  et  al.,  1970);  ZB  is, technically,  the  ratio  at  the  Earth's  surface 
between  the  horizontal  component  of  an  oscillating  electric  Held  and  the  horizontal  transverse 
component  of  the  magnetic  field  that  appears  in  response  to  it,  and  on  grounds  of  symmetry  we 
have  in  this  case  that  this  total  magnetic  induction  Bf  must  be  perpendicular  to  the  induced 
electric  field  £q&  (as  is  usually  the  case  in  general),  yielding  the  field  direction  indicated  by 
(2.11)  (arguments  of  scaling  and  symmetry  suggest  that  if  Sr  has  a  vertical  component  then  it 
must  be  much  smaller  in  magnitude  than  the  horizontal  component).  At  a  frequency  of  one 
cycle  per  hour,  in  a  typical  measurement  Z£  might  have  a  magnitude  of  0.00025  ft,  and  a  phase 
of  about  45*.  Br  is  the  sum  of  Bf  and  Bp,  the  slowly  fluctuating  magnetic  induction  field  that 
originally  induced  B},  and  by  the  previous  argument  relating  to  the  magnetic  induction  number 
we  have  that  6jr  must  be  essentially  uniform  in  space  throughout  the  ocean  and  f  or  some  dis¬ 
tance  above  the  ocean  surface;  accordingly,  we  get  Bj-  over  this  altitude  range  simply  by  adding 
to  Bf  whatever  is  necessary  to  bring  its  value  at  the  ocean  floor  up  to  that  given  by  (2.10). 
Thus,  S/r  comes  to  . 

j  ,  (2.11) 

and  the  expression  for  Br  is 

^(z0)  -  mo£o.P  |  -  K*o>  J  *  (212) 

where 


<r(z)dz  ,  0<z0<£ 


(2.13) 


and 

X(r0)-ZU»  ,  z0>D 


(2.14) 


For  an  ocean  4  kilometers  deep,  for  an  inducing  field  oscillating  at  one  cycle  per  hour,  t(D)  is 
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smaller  than  \JZ£  by  a  factor  of  about  3. 


fx 
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3.  The  perterblng  effect  of  shallow  waves  at  the  surface 

We  now  suppose  the  surface  to  be  described  by  the  expression 

l- a  exp(ilcH-X)  e-,ul  ,  (3.1) 

where  (  is  the  height  of  the  surface  above  its  mean  level,  a  is  the  wave  amplitude,  kH  is  the 
horizontal  component  of  the  wavevectoi ,  which  for  convenience  we  will  express  in  the  form 

—  i It  (cos6,  sin0, 0)  ,  (3.2) 

where  k  is  real  and  positive,  H  is  the  position  vector,  3?«(v,y*c),  and  <o  is  the  angular  fre¬ 
quency  of  the  wave.  In  practical  applications  the  surface  is  considered  to  be  given  by  either  the 
real  or  imaginary  part  of  expression  (3.1),  whichever  is  more  convenient,  and  the  resultant 
magnetic  field  is  determined  by  taking  the  same  part  of  the  complex  field  expression  derived 
from  (3.1). 

In  order  to  meet  the  boundary  condition  that  7  at  the  surface  should  have  no  com¬ 
ponent  normal  to  the  surface,  it  is  necessary  to  add  a  perturbation  term  to  the  uniform  current 
distribution  initially  assumed  to  exist  beneath  the  surface.  Furthermore,  since  the.,  current  dis¬ 
tribution  within  the  ocean  is  ultimately  induced  by  whatever  electric  field  is  present, 

7«or  2?  ,  (3.3) 

where  1  is  the  electric  field  and  <r  is  the  electrical  conductivity  of  seawater  (which  may  be  con¬ 
sidered  for  our  purposes  to  be  a  constant,  as  the  depth  range  within  which  the  surface  waves 
exert  a  noticeable  electrical  influence  is  too  thin  to  exhibit  a  significant  conductivity  variation), 
7  may  be  viewed  as  the  gradient  of  a  scalar  function; 

7  -V*  .  (3.4) 

If  we  make  the  reasonable  assumption  that  separation  of  electrical  charge  within  the  body  of  the 
ocean  is  negligible  for  all  m  of  interest,  we  also  have 

V’7-0  ,  (3.5) 

which  with  (3.4)  implies 

VJ*  -  0  .  (3.6) 

The  only  solutions  to  (3.6)  that  can  conveniently  be  reconciled  with  (3.1)  have  the 

form 
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<f>  -  (  A  +  B  expiiicxll) }  e~lml  ,  (3.7) 

where 

Tc\  -  k  (cosfi,  sin5,~/)  (3.8) 

tod 

£2-  *  (aw#,  sirtf,  /)  ;  (3.9) 

the  first  term  of  (3.7)  dor.ps  exponentially  downward,  and  the  second  term  grows  exponentially 
downward.  For  an  ocean  of  finite  and  uniform  depth  the  coefficients  A  and  B  can  be  related  by 
use  of  the  boundary  condition  at  the  ocean  floor  that  the  electrical  current  flow  at  this  interface 
has  no  vertical  component,  but  as  long  as  the  order  of  magnitude  of  the  wavelength  of  the  sur¬ 
face  wave  is  less  than  that  of  the  depth  of  the  ocean,  which  for  our  purposes  will  generally  be 
the  case,  the  second  term  of  (3.7)  is  negligible  compared  to  the  first;  accordingly,  to  simplify 
the  calculations  we  will  henceforth  assume  that  B-0  and  take  giving  a  current  distribu¬ 
tion  beneath  the  surface  of 

7**  JoSt  +  lc  IA  exp iM'Tt)  e~lmt  ,  (3.10) 

£  -  it  (coafl,  sinfl, -/)  .  (3.11) 


The  value  of  A  may  be  determined  from  the  boundary  condition  at  the  surface  that 

7-7-0  atz-l  ,  (3.12) 

where  7  is  a  vector  normal  to  the  surface;  such  a  vector  is 


7- 


_«i.  _!£.  1 

Bx'  By' 


(3.13) 


If  (3.10)  is  substituted  into  (3.12)  and  terms  of  second  and  higher  orders  in  Cka)  are  discarded, 
one  gets 


A  -  k/ocosfl  (3.14) 

under  the  condition  that 

\ka  |  «  1  .  (3.15) 

This  gives  from  (3.10) 

7 -  J0  |  Jc  -  a  cost)  exp  (ifc-J)  e f  )  .  (3.16) 
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hi  general,  the  magnetic  held  B  OTo)  due  to  an  electrical  current  distribution  within  a 
volume  V  is  given  by 


(3.17) 


(Reitz  and  Milford,  1967,  p.154);  for  this  problem  this  integral  takes  the  form 

■•'-si'W-T?  • 


(3.18) 


where  7 Of)  is  os  given  by  (3.16).  It  is  convenient  for  purposes  of  calculation  to  express  B  as 
the  sum  of  B<{  and  By , 


B  (?«,)  -  3i  OT0)  +  3V  0?o)  , 


(3.19) 


where  is  the  field  due  to  the  part  of  the  current  distribution  that  is  between  z*--D  and  z~0 
and  Bs  is  the  field  due  to  the  part  of  the  current  distribution  between  z— 0  and  z-£,  and  to 
break  'By  down  into  Byv  and  By? , 


By  (?q)  •*  Byy  C?0)  +  Byp(Xo)  ,  (3.20) 

where  is  due  to  the  uniform  component  of  the  current  distribution,  JqSc,  and  Byp  is  due  to 
the  perturbation  thereon  which  constitutes  remainder  of  the  distribution.  Bvu  has  already  been 
determined  in  section  2. 

Byy  is  given  by 

c  "  °- 


SrnOfo) 

-as  -■  h» 

(cos0,  sin0,  -/)  x  0?-3?o) 


(3.21) 


C  —  no  Jo  a  k  cosfl  , 


exp  (/fo?) 


(3.22) 


where  we  have  simplified  the  integration  by  extending  the  lower  limit  of  integration  from  z--D 
down  to  z—-oo;  provided  that  the  wavelength  of  the  surface  wave  is  somewhat  less  than  the 
depth  of  the  ocean,  as  has  already  been  assumed,  the  error  introduced  by  this  extension  should 
be  negligible.  Making  the  substitutions 


u  «*  (x-xo)  cos0  +  (y-yo)  sind 
v  ■*  (y-yo)  cos#  -  (x-xq)  sin# 


(3.23) 

(3.24) 
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w  -  -  (x-Zq) 

Uq  —  *o  cos0  +  yo  sin® 
tt  —  x  cos 0  +  £  sine 
5  -  .P  cose  -  x  sine 
turns  (3.21)  into 

ctto °e**°  7  .  7 


(3.25) 

(3.26) 

(3.27) 

(3.28) 


5np  C?g)  * 


4ir 


jf  </u  /  dv 


(3.29) 


*o 


(/8+2)v+  P(w-/u)  -** 

(ui+vW)*2  tf  C 


where  x0  is  the  altitude  of  the  observation  point  relative  to  the  ocean  surface.  The  term  in 
( ffi+2 )  does  not  contribute  to  the  integral,  as  it  gives  an  integrand  antisymmetric  in  v.  Further 
reduction  of  the  integral  is  aided  by  the  identity 


f  <fr  2 

(x^fc)*'2  e 


(3.30) 


(Orobner,  sec.  213,  eq.  1,  p.  34);  use  of  this  identity  in  the  integration  over  v  gives  frum 
(3.29) 


-M  In 


gtku  gkw 


*0 


w+iu 


(3.31) 


which  yields  easily  to  contour  integration  in  u  (over  the  upper  half  plane)  followed  by  standard 
integration  in  w,  to  give 


EypOto)  -  V4  MO  h  a  C0S(9  rlto°'  e“'«' 
x  (-sing,  cos0,  0)  . 


(3.32) 


In  the  computation  of  5s  the  perUnbation  part  of  the  current  distribution  contributes 
only  in  second  order  in  (ka),  and  thus  can  be  ignored,  leaving  only  the  uniform  component  of 
the  current  distribution  to  be  considered.  The  appropriate  integral  expression  is  then 


(3.33) 
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or 


Ss  <X(y)  -  .  f  dx  f  dy  \dz 

™  i  i,  o 

(3.34) 

..  {0,  (T-Zoi.-O’-^o)} 

\X-*o\3 

The  small  range  of  integration  of  z  in  this  integral  suggests  a  Taylor  expansion  of  the  integrand 
over  this  variable;  if  this  is  done,  only  the  first  ,  term  (the  constant  term)  contributes  to  in 

first  order  in  ka.  Doing  this,  we  get 

tsw-zg-l* 

(3.35) 

10.  *o.  (v->o»  .  , 

Using  the  expression  (3.1)  for  (  and  making  the  variable  changes  (3.23),  (3.24),  and  (3.26) 

gives  from  this 

,,  J  a  J*0  f-lmt  "L  “ 

isW— Jdu  Jd, 

(3.36) 

(0,  z0, u  sin9+v  cos0)  lku 
*  (u^fv^)*1  *  ’ 

• 

from  which  the  integration  over  v  gives 

"rt  /„  \  no  Jo  a  e~imt  f  (0,  z0,  u  sin0)  e 

^  2rr  J.*  (u+tzoHu-izo)  ' 

(3.37) 

which  yields  to  contour  integration  in  u,  closing  the  contour  in  the  upper  half  plane,  to  give 

<20) - 'h  Ho  Jo  o  «ta°  e“lteo'  rM  T (zo) 

(3.38) 

J  (0, 1,  /  sinfii)  ,  zo>  0 
(0,-1,/  sing)  ,  z0<  0  . 

(3.39) 

Combining  this  with  (3.31)  gives  for  2? w,  the  field  due  to  the  perturbing  effects  of  the 
surface  wave, 

5*  (5?o) "  —  VS  Mo  o  elku°  e-|teo1  e~lml  T(zq) 


(3.40) 


t 
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t 


i 

1 


l 

i 

\ 

i 

t 

« 

t 

i 

i 

i 

t 

r 


7<2o) 


sind  (cosd,  sinfl,  0  ►  *o  >  0 

sind  <«wd,  sind,  /)  -  (0,  2,  0)  ,  *o  <  0  , 


(3,41) 


to  which  the  field  computed  in  section  2  for  a  uniform  current  distribution  within  a  flat-surfaced 
ocean  must  be  added  to  give  the  complete  magnetic  held  to  first  order  in  (ka).  As  a  check, 
(3.40)  and  (3.41)  have  been  verified  by  use  of  the  curl  expression  (2.S).  Note,  as  a  special 
case,  that  when  the  wave  crests  are  perpendicular  to  the  x  direction  there  is  no  wave-induced 
contribution  to  the  field  above  the  surface.  This  is  reasonable  on  physical  grounds,  as  in  this 
case  the  total  vertically  integrated  horizontal  electric  current  flow  within  the  ocean  is  unchanged 
from  the  flat  surface  condition  and  unchanging  with  position  everywhere  in  the  ocean;  there¬ 
fore,  by  Ampere's  law,  the  difference  AB  between  the  magnitudes  of  the  horizontal  component  , 
of  the  magnetic  field  above  the  ocean  surface  and  beneath  the  ocean  floor  is  unchanged  from 
the  flat  ocean  case,  and  considerations  of  symmetry  rule  out  a  vertical  field  component  for  this 
geometry.  By  contrast,  when  the  wavr  crests  are  parallel  to  £  there  is  more  total  electrical 
current  flow  through  the  water  beneath  a  crest  than  beneath  an  adjacent  trough,  and  hence  AB 
is  variable  from  place  to  place. 


4.  Application  to  Adds  induced  by  bottom  topography 

The  work  in  the  previous  section  can  be  applied,  with  minor  modifications,  to  the  cal¬ 
culation  of  magnetic  anomalies  produced  by  certain  types  of  bottom  features.  Any  reasonably 
smooth  bottom  feature  can  be  described  by  a  series  of  terms  of  the  form  of  (3.1)  with  a>  set  to 
zero,  and  the  problem  is  then  reduced  to  that  of  the  previous  section  turned  upside-down.  The 
usefulness  of  this  approach  is  limited  in  that  the  lateral  dimension  of  the  feature  to  be  investi¬ 
gated  must  be  somewhat  greater  than  its  vertical  dimension  in  order  to  meet  the  small-scale 
approximation  (3.15)  used  in  developing  the  mathematics.  This  application  is  discussed  in 
detail  in  section  6  of  chapter  2. 

In  the  event  tha  .e  lateral  dimension  of  the  feature  is  on  the  order  of  the  depth  of 
the  ocean,  both  terms  of  the  general  solution  (3.7)  for  the  electrical  current  potential,  that 
which  increases  exponentially. with  depth  and  that  which  decreases  exponentially  with  depth, 
must  be  used,  and  the  subsequent  steps  of  the  derivation  altered  to  suit.  Since  the  vertical 
scale  of  the  bottom  feature  wilt  generally  be  much  greater  than  the  amplitudes  of  whatever 
waves  there  are  on  the  ocean  surface,  it  is  a  reasonable  simplifying  assumption  in  this  develop¬ 
ment  to  consider  the  ocean  surface  to  be  completely  fiat. 
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5.  Surface  wares  moving  through  the  Earth’s  magnetic  held  f 

Surface  waves  will  be  accompanied  by  motion  of  the  seawater  beneath  the  surface,  and  i 

the  movement  of  charge-carrying  fluid  through  the  Earth's  magnetic  field  will  result  in  an 
induced  electrical  current  distribution;  this  current  distribution  will  in  turn  give  rise  to  a  mag¬ 
netic  field  modifying  the  Earth’s  field.  This  subject  has  been  covered  by  Crews  and  Futterman 
(1962),  Warburton  and  Caminiti  (1964),  and  Weaver  (1965),  but  for  the  sake  of  completeness 
an  abbreviated  derivation  of  the  effect  is  given  here  in  the  mathematical  framework  developed  f 

in  sections  2  and  3  of  , this  chapter.  | 


The  current  distribution  induced  by  the  water  movement  is  given  by 
7“<rTfx?  +  7c  , 


(5.1)  • 


where  a  is  the  electrical  conductivity  of  the  water,  taken  here  to  be  constant,  7  is  the  water 
velocity  distribution,  and  T  is  the  Earth's  magnetic  field,  which  may  be  taken  to  be  locally  uni¬ 
form.  7c  is  *  correction  term  whose  purpose  is  to  insure  that  the  boundary  condition  is  met  that 
there  be  no  electrical  current  flow  through  the  surface.  If  we  assume  the  water  flow  to  be  irro- 
tational  (for  a  discussion  of  the  validity  of  this  assumption  see  Lamb,  1945,  sec.  33,  pp.  35-37), 
then  the  condition 


v-7-o 

is  automatically  satisfied  provided  that  the  same  condition  holds  for  7c  ; 
V-7-o-  V'(Vx?)  +  V-7C 

-  a  (^,•CVx7>-V•<Vx?,)}  +  V-7C 

-0  . 


(5.2) 


(5.3) 


Another  consequence  of  the  assumption  of  irrotational  water  flow  is  that  the  velocity  field  may 
be  viewed  as  the  gradient  of  a  scaler  function,  and  if  we  also  suppose  that  the  seawater  is 
incompressible,  which  is  safe  enough,  we  have 


(5.4) 


V’V  «*•  0  , 

If  4>  is  the  scalar  function,  then 

T-V*  (5.5) 

Vty-0  .  (5.6) 

The  solution  of  (5.6)  that  most  conveniently  fits  with  the  surface  equation  (3.1)  has  the  form 


iii 
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#* 

1 


apiUc-n)  *“**  , 

where  H  is  as  given  in  (3.1 1),  and  if  we  make  the  constraint  that 

an  approximation  that  is  correct  to  first  order  in  (ka),  we  find  that 


"--'If  • 


V  -  (ensg,  sinS,  —i)  (a  a  expUJcll)  t ~t*‘l  . 


7  is  then  given  by 


7— <r  w  a  expiry)  «'~M  (cosS,  sing,  -/)  x  T 
+  7c  . 


(5.10) 


(5.11) 


7C  will  have  zero  divergence  as  required  if  it  is  proportional  to  V,  and  the  boundary  condition 
that  there  be  no  electrical  current  flow  through  the  surface,  given  by  (3.12),  is  satisfied  to  first 
order  in  (ka)  if  7,-0;  the  result  is 


7cmm-l  (cosfi,  sinS.-O  cr  <o  a  ( («w9,  sin®,-/)  x  ?}, 
« 

x  cg>  dEx)  *"lml 

7—  <r  6i  a  x  bb*0,  --cos 9,  0)  expiiV'X) 

X  -  <  (coitf,  sinS,-/)-7  . 


(5.12) 


(5.13) 

(5.14) 


Ihe  magnetic  field  arising  from  this  current  distribution  may  be  found  by  inserting 
(5.13)  into  (3.18).  As  in  section  3,  we  find  it  convenient  to  divide  the  integral  into  an  integral 
over  the  volume  below  z-0  and  an  integral  covering  the  complementary  region  between  z-0 
and  the  surface  at  upon  doing  this,  we  find  by  inspection  that  the  second  integral  contri¬ 
butes  nothing  below  the  second  order  in  (ka)  to  the  magnetic  field,  leaving  only  the  first 
integral  to  evaluate.  Making  the  substitutions  (3.23)  through  (3.26)  and  simplifying  the 
integration  by  assuming  infinite  depth  for  the  ocean  as  in  section  3  gives  the  integral 


50?<,)  -  -•£-  Jciu  Jdv  Jdw 

om  -m  r0 

(w>  cosfi,  w  sine,  u)  lku  kw 

<«2+vW)y2  *  * 


(5.15) 


if**- 
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C  «■  fio  <r  tt  a  x  t1*9  efa° 


Integration  over  v  gives  from  this 

BW-  -£1*  ]*»  (w  •»  r»  , 


(5.16) 


(5.17) 


from  which  contour  integration  in  u  (closing  the  contour  in  the  upper  half  plane)  followed  by 
integration  over  w  gives 


5(*o)  IMmoo’  a  X  t>h>0  ?(«o) 


t(r  o)- 


<COS0,  sin®,  /) 


“|  (coa9,  sin 9,1)  +  2uzq  ( cosO ,  sirtf,  -/) 


,  *o>  0 
,  *o  <  0  • 


This  has  been  verified  with  the  curl  expression  (2.5). 


(5.18) 


(5.19) 
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6.  Gradients  of  the  held  components 

Above  the  surface,  the  field  derived  in  section  2  for  a  uniform  electric  current  flow, 
given  by  (2.2)  through  (2.4),  is  completely  uniform  itself,  and  therefore  has  a  zero  gradient. 
Below  water,  though,  the  field  magnitude  decreases  linearly  with  increasing  depth;  the  gradient 
of  the  y-component  of  the  field,  the  only  nonzero  component,  is  given  by 

of  .  (6.1) 

The  calculation  of  the  tensor  gradients  of  the  magnetic  field  derived  in  section  3,  and 
given  by  (3.37)  and  (3.38),  is  made  easy  by  the  fact  that  all  of  the  field  components  have  the 
same  spatial  dependence, 

/0?) -«*<£*)  ,  (6.2) 

where  Tc  is  given  by  (3.8)  where  z<0  and  by  (3.9)  where  z>0;  the  gradient  of  this  expression 
is  given  by 

V/CP)-/E/C?)  ,  (6.3) 

where  Tc  is  defined  as  above.  Accordingly,  the  appropriate  component  gradients  may  be  derived 
just  by  multiplying  the  individual  field  components  by  He  with  Tc  given  by 

|  k  (cos9t  sit#,  /)  ,  z  >  0 

*  "  J  *  (cos9,  sin 0,-1)  ,  *  <  0  .  • (6l4) 

Examining  (5.18)  and  (5.19),  the  field  expression  derived  in  section  5,  we  seo  that  the 
field  above  the  surface  has  the  spatial  dependence  (6.2),  and  the  gradients  of  its  vector  com¬ 
ponents  are  found  in  the  same  way  as  those  for  the  field  derived  in  section  3.  Beneath  the  sur¬ 
face  the  field  expression  has  two  terms;  one  of  these  has  the  spatial  dependence  (6.2),  but  the 
other  has  the  spatial  dependence 

g(X)-z  expUTcX)  ,  (6.5) 

vhoue  gradient  is 

V*0?)-  ff+i  ,00  ,  (6.6) 

7  I 

where  V  is  defined  as  in  (6.4)  for  z<0.  In  order  to  calculate  the  gradient  of  a  given  underwater 
field  component,  one  multiplies  the  appropriate  component  of  the  first  term  of  the  field  expres¬ 
sion  by  He ,  then  multiplies  the  corresponding  component  in  the  second  term  by  (M-’+z/z),  and 


^  ■  !•  M  «  «  ►  ■  *M  ,*\  <•  .  *  k  *.  -  »  k 
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finally  adds  together  the  resulting  vectors. 


18 


7.  Discussion 

In  the  open  ocean,  the  ionospherically  induced  electric  fields  that  one  might  expect  to 
find  on  an  electromagnetically  quiet  day  will  have  fluctuations  on  the  order  of  1  p.  V/m,  induced 
mainly  by  the  solar  daily  variation,  whose  amplitude  will  be  roughly  2x  1CT8  tesla,  whereas  dur¬ 
ing  an  Intense  magnetic  storm  the  field  intensity  might  go  as  high  as  10  n  V/m ;  in  the  neigh¬ 
borhood  of  a  coastline,  where  there  are  boundary  effects,  the  field  intensities  will  typically  be 
larger  than  those  of  the  open  ocean  by  a  factor  of  2  to  3.  Assuming  a  uniform  electrical  conduc¬ 
tivity  for  the  ocean  of  3.3  S/m,  this  implies  induced  electrical  current  densities  in  the  open 
ocean  of  from  3.3  fiA/m2  up  to  33  ixA/m J,  with  larger  values  near  a  coastline.  Use  of  this 
range  of  current  density  values  in  equation.  (2.5),  assuming  an  ocean  depth  of  5.0  kilometers, 
gives  for  Bi,  the  magnetic  induction  in  the  air  over  the  ocean  due  to  induced  electrical  current 
flow  within  it,  a  magnitude  ranging  from  l.Ox  10“'  tesla  to  l.Ox  Iff*7  tesla  and  a  direction  hor¬ 
izontally  perpendicular  to  the  subsurface  electrical  current  flow.  By  comparison,  the  Earth's 
magnetic  field  over,  the  open  ocean  off  of  La  Jolla,  California,  is  roughly  6x  10-5  tesla  in  inten¬ 
sity  and  has  a  tilt  of  about  60s  from  the  horizontal. 

The  magnitude  of  the  magnetic  induction  generuted  by  the  interaction  of  an  ocean 
wave  with  the  uniform  current  density  varies  with  the  angle  that  the  wave  crest  makes  with  the 
uniform  Bow,  but  a  root  mean  square  average  of  (3.39)  over  time,  lateral  position,  and  this 
angle  gives  for  z>0 

|5|  —  ftftoJoa  exp(r-2irz/\)  ,  (7.1) 

where  a  is  the  wave’s  amplitude  and  X  is  its  wavelength.  A  typical  wave  amplitude  on  a  calm 
day  is  roughly  1  meter,  and  during  a  storm  one  might  see  a  typical  wave  amplitude  on  the  order 
of  5  meters.  Assuming  a  wave  amplitude  of  1.0  meter  and  using  the  current  density  range  of 
from  3.3  pA/m1  to  33  fiA/m1  gives  at  just  above  the  surface  a  magnetic  induction  magnitude 
range  of  from  2.1X10-12  tesla  to  2.1x10-"  tesla,  assuming  that  the  wave  crests  are  parallel  to 
the  electric  current  flow. 

Performing  the  same  averaging  process  on  (5.18),  we  gel  for  the  mean  magnitude  in 
the  air  of  the  magnetic  induction  of  surface  waves  moving  through  the  Earth’s  field 

|5|  -  'h^tr  a  (talk)  G  exp(-2irz/\)  (7.2) 

G  -  V  (tf+Vitf)  ,  (7.3) 

where  F,  is  the  vertical  component  of  the  Earth’s  field  and  Fh  is  the  magnitude  of  the  horizon¬ 
tal  component.  Most  of  the  surface  wave  energy  in  the  open  ocean  off  of  La  Jolla  falls  in  the 


; 
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range  of  wave  periods,  from  5  seconds  to  IS  seconds,  equivalent  to  an  angular  frequency  range 
of  from  0.42  .ate-1  to  1.26  arc"1;  using  the  dispersion  relation  for  surface  waves  on  a  deep 
ocean, 

^wmgic  (7.4) 

(SommerfelcL,  pp.  168-173),  where  g  is  the  gravitational  acceleration  and  has  a  value  of  about 
9.81  m/sec *,  we  calculate  a  range  of  phase  velocities  (u/k)  of  from  7.8  m/sec  to  23.4  m/sec. 
The  value  for  O  in  the  vicinity  of  La  Jolla  comes  to  about  6x  10" 5  tesla,  and  if  we  again  assume 
a  wave  amplitude  of  1.0  meter  and  an  electrical  conductivity  for  the  ocean  of  3.3  S/m  we  get  a 
range  of  magnitudes  for  the  magnetic  induction  at  just  above  the  ocean  surface  of  from 
6x  10“*°  tesla  to  2x  10"9  tesla. 

The  uniform  held  discussed  in  the  first  paragraph,  that  is  produced  by  the  ionospheri- 
tally  induced  current  flow,  is  at  the  weaker  calculated  extreme  a  little  less  than  four  orders  of 
magnitude  down  from  the  Earth’s  field,  and  will  generally  not  be  detectable.  The  two  wave 
interaction  effects  are  roughly  five  orders  of  magnitude  down  from  the  Earth's  field,  but  their 
periodic  nature,  both  spatially  and  temporally,  should  allow  them  also  to  be  detected.  Except  in  • 
the  presence  of  e  powerful  magnetic  storm  the  wave-magnetic  field  interaction  effect  discussed 
in  the  previous  paragraph  will  typically  be  roughly  two  orders  of  magnitude  stronger  than  the 
wave-current*  interaction  effect  discussed  in  the  second  paragraph,  making  the  latter  difficult  to 
resolve  from  the  former,  as  above  the  surface  they  have  the  same  spatial  dependence. 
Although  these  effects  both  scale  linearly  with  the  wave  amplitude,  making  their  ratio  indepen¬ 
dent  of  this  quantity,  the  wave-magnetic  field  interaction  effect  is  also  directly  proportional  to 
the  phase  velocity  of  the  wave  in  question,  which  increases  linearly  with  the  wave’s  period  as 
long  as  the  deep  water  approximation  is  valid.  The  two  wave  interaction  effects  are  seen  from 
rough  calculations  to  have  comparable  magnitudes  in  the  extreme  case  of  a  surface  wave  with  a 
period  of  (U6  second  or  less  in  the  presence  of  a  strong  magnetic  storm. 


Chapter  II 


Magnetic  Anomalies  Induced  by  Bottom  Features 


1.  Introduction 

We  consider  in  this  chapter  magnetic  anomalies  associated  with  the  presence  of  rocky 
projections  from  and  depressions  in  an  otherwise  level  ocean  floor.  Only  magnetic  effects  due  to 
electrical  current  flow  within  the  water  of  the  ocean  are  considered;  magnetic  effects  due  to 
magnetization  of  the  ocean  floor,  to  the  magnetic  properties  of  large  projections  of  magnetic 
minerals  from  the  ocean  floor,  and  to  irregularities  of  magnetization  of  the  seafloor  rocks  are 
covered  in  other  works  (see  for  example  Larson  et  al.,  1974,  or  Vacquler,  1972).  In  this  treat¬ 
ment,  a  large  scale,  horizontal  uniforni  electric  current  flow  driven  by  sources  of  EMF  outside 
of  the  ocean  is  assumed,  and  the  local  distortion  of  this  flow  by  the  presence  of  a  bottom 
feature  (trnted  as  an  insulator),  is  determined;  then,  by  means  discussed  in  section  2,  the 
magnetic  field  perturbation  associated  with  this  electrical  current  flow  distortion  is  calculated. 
As  the  scale  of  the  bottom  feature  in  question  will  typically  be  several  orders  of  magnitude 
greater  than  the  amplitudes  of  any  waves  on  the  ocean  surface,  the  simplifying  assumption  of  a 
completely  flat  ocean  surface  will  be  made. 

In  section  s,  in  order  to  give  some  feeling  for  orders  of  magnitude,  numerical  calcula¬ 
tions  arc  made  of  magnetic  fields  and  the  tensor  gradients  thereof  associated  with  a  submerged 
hemisphere  with  a  radius  half  the  depth  of  the  ocean  and  resting  on  a  level  ocean  floor.  In  sec¬ 
tion  4,  the  problem  of  a  localized  bottom  feature  of  arbitrary  shape  is  considered,  and  section  5 
gives  a  quick  algorithm  for  approximating  the  electrical  current  flow  deflection  by  bottom 
features  of  small  aspect.  Section  6  covers  in  detail  the  magnetic  effects  of  small  features  of  low 
aspect  on  the  ocean  floor  (such  as  a  low  ridge  or  shallow  trough),  a  subject  touched  lightly  upon 
in  section  4  of  chapter  \.  Finally,  section  7  briefly  discusses  the  problem  of  the  magnetic 
effects  due  to  the  interaction  of  seawater  current  flow  over  a  bottom  feature  with  the  Earth’s 
magnetic  field. 
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2.  Mathematical  methods 


2.1  The  basic  Integration  formula 

The  basic  formula  for  computing  the  magnetic  field  3  due  to  an  electrical  current  den¬ 
sity  distribution  7  within  the  ocean  is  the  law  of  Biot  and  Savart  as  applied  to  a  continuous 
electrical  current  density  distribution, 


(2.1) 


7(x)  is  assumed  to  have  already  been  adequately  determined  for  practical  purposes,  and  may  be 

taken  as  a  given. 

In  practical  calculations,  the  selected  volume  of  integration  is  broken  down  into  rec¬ 
tangular  boxes,  and  the  electrical  current  density  vector  7(xf7  and  its  tensor  gradient  at  the 
geometrical  center  jg  of  each  box  are  determined;  the  current  density  distribution  throughout 
the  box  is  then  assumed  to  be  given  by  th?  formula 


Jj<*) -/,(£)+  0r-$)V/,05) 


(2.2) 


Also,  it  is  assumed  that  the  distance  between  x0  and  £  is  much  preator  than  the  maximum 
dimension  of  the  box,  allowing  the  expansion  of  the  denominator  of  the  integrand  of  (2.1)  in  a 
binomial  series.  Oiven  these  assumptions,  one  can  write  an  algebraic  expression  for  the  contri¬ 
bution  to  the  integral  (2.1)  and  its  tensor  gradient  of  any  given  box  in  terms  of  the  box  dimen¬ 
sions,  the  position  J£  of  the  center  of  the  box,  and  the  electrical  current  density  and  its  tensor 
gradient  at 

We  start  with  a  variation  of  (2  1), 
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where 
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ro  -  Xo  -  £  , 


(2,4) 


?  ,  (2.5) 

and  R  denotes  the  rectangular  volume  of  integration  centered  on  the  position the  vector  .v! 
is  a  small  ofTset  from  xjj  in  an  arbitrary  direction,  and  is  used  to  find  the  tensor  gradient  of  £  at 
Xq.  For  the  current  density  distribution  in  the  numerator  of  the  integrand  we  write 
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and  expansion  of  the  denominator  gives 
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Insertion  of  (2.6)  and  (2.8)  into  (2.3)  gives  a  series  expansion  in  'rj,  xj,  and  V  for  the 
integrand,  only  a  few  terms  of  which  arc  of  practical  interest  to  us.  The  terms  in  this  series 
that  are  linear  in  xl  are  sufficient  for  the  determination  of  the  tensor  gradient  of  the  magnetic 
Held,  and  so  terms  quadratic  or  higher  in  x\  may  be  ignored;  also,  due  to  the  symmetry  of  the 
volume  of  integration  about  j£,  the  integrals  of  terms  linear  or  cubic  in  V  vanish,  and  those  of 
terms  quartic  or  of  higher, powers  in  V  will  in  general  be  small  enough  in  comparison  to  the 
integrals  of  terms  quadratic  or  constant  in  V  to  be  safely  disregarded  in  practical  calculations. 
These  two  conditions  together  eliminate  all  terms  in  the  expansion  (2.8)  derived  from  terms 
quartic  or  higher  in  z,  and  most  of  the  terms  in  smaller  powers  of  z  in  this  expansion  When 
all  terms  of  interest  are  taken  into  account,  we  have  from  (2.3) 
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where 


(2.21) 

(2.22) 


(2.23) 


dt  being  the  dimension  of  the  volume  of  integration  in  the  direction  Jr, .  The  basic  relation  used 
for  the  evaluation  of  the  above  integrals  is 

J^vtvjdi r  -  rjVa  8r$  ,  (2.24) 


where  6y  is  the  JCronecker  delta.  Expression  (2.1 1)  may  be  rewritten  in  the  form  , 

2f(xetxj>-  B(xJ)  +  KJ  ,  (2.25) 

where  the  construction  of  the  tensor  7  is  just  a  matter  of  sorting  the  coefficients  of  the  vector 
elements  of  jcJ  in  (2.11)  into  the  proper  tensor  elements  of  T.  When  the  field  S  is  expressed 
in  this  form,  it  is  apparent  that  T  is  the  tensor  gradient  of  S  at  the  point  x^. 


2.2  CorrectUo  for  edge  effect  errors 

One  difficulty  in  the  practiced  use  of  (2.1)  is  that  it  makes  the  implicit  assumption  that 
the  entire  electrical  current  density  distribution  is  contained  within  the  volume  of  integration, 
and  use  of  a  volume  of  integration  that  has  a  current  flow  through  any  part  of  its  surface  results 
in  the  appearance  of  nonphysical  elements  in  the  results  of  the  integration.  As  an  example  of 
this  effect  mo  consider  the  case  of  a  semi-infinite  wire  extending  from  the  origin  in  the  positive 
St  direction  out  to  infinity  and  carrying  an  electrical  current  in  that  direction,  in  this  case  (2.1) 
becomes 
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this  reduces  to 
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where 


(Cxr**+P«s! 


(2.31) 
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(CRC,  p.  322,  eq.  196).  For  large  distances  from  the  origin  and  relatively  much  smaller  dis¬ 
tances  from  the  origin  this  gives 
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the  field  around  an  infinitely  long  wire,  as  one  would  expect  on  physical  grounds.  However,  if 
we  use  Ampere's  law  to  compute  from  (2,29)  the  electrical  current  density  outside  of  the  wire, 
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which  on  physical  grounds  should  be  zero,  we  get 

70?)  -  —  VxSC?)  (2.34) 
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the  physical  interpretation  of  this  expression  indicates  a  current  density  distribution  converging 
on  the  end  of  the  wire  at  the  origin  in  such  a  way  as  to  feed  a  total  current  /o  into  the  wire  end. 
In  general,  it  seems  that  assuming  a  source  or  sink  of  electrical  current  in  the  postulated  electri¬ 
cal  current  density  distribution  results  In  a  mathematical  artifact  that,  if  taken  at  face  value, 
cancels  out  the  source  or  fills  in  the  sink. 

In  the  calculation  of  magnetic  fields  associated  with  localized  bottom  features,  one  will 
typically  be  dealing  with  an  electrical  current  density  distribution  anomaly  the  bulk  of  whose 
magnetic  effects  are  produced  by  the  part  of  the  distribution  in  the  immediate  neighborhood  of 
the  bottom  feature,  but  which  gives  minor  contributions  to  the  magnetic  field  at  substantial  dis¬ 
tances  from  the  feature;  in  numerical  calculations  of  the  magnetic  field  the  approach  that  first 
suggests  itself  is  to  consider  only  the  volume  immediately  about  the  bottom  feature,  tolerating 
errors  of  perhaps  a  few  percent  in  the  field  values  of  Interest  in  order  to  avoid  expensive 
integration  over  large  volumes  of  seawater.  On  the  basis  of  the  example  of  the  semi-infinite 
wire,  it  appears  that  if  this  approach  is  used  it  must  be  used  with  care,  as  neglecting  the 
extended  part  of  the  current  density  distribution  not  only  results  in  errors  of  omission,  blit  may 
also  contaminate  the  calculated  values  with  mathematical  artifacts. 

A  method  for  dealing  with  the  artifact  problem  may  be  derived  by  viewing  an  infinitely 
long,  straight,  current  carrying  wire,  for  which  there  is  an  artifact-free  solution  for  the  magnetic 
field  and  its  tensor  gradients,  as  being  composed  of  two  semi-infinite  wires  pointing  outward 
from  the  origin  in  opposite  directions  and  carrying  currents  of  equal  magnitude  and  opposite 
sign.  From  the  observation  that  the  artifacts  associated  with  the  two  wire  segments  must  exactly 
cancel,  we  deduce  that  any  terms  in  the  magnetic  field  expression  for  a  semi-infinite  wire  that 
contribute  to  the  artifact  must  be  asymmetric  with  respect  to  simultaneous  reversal  of  wire 
direction  and  current  direction,  and  we  infer  from  symmetry  arguments  that  any  singularities  in 
the  artifact  terms  must  be  at  the  origin.  Given  that  these  constraints  are  sufficient  to  uniquely 
identify  those  terms  contributing  to  the  artifact,  once  an  artifact  correction  expression  for  a 
semi-infinite  wire  has  been  derived  it  may  be  use  as  a  Green’s  function  for  a  general  current 
density  flow  through  a  surface. . 
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To  derive  an  expression  Tor  a  semi-infinite  current  carrying  wire  pointing  from  the  ori¬ 
gin  in  an  arbitrary  direction,  we  start  with  a  binomial  expansion  of  t2,l),  after  the  manner  of 
(2.3); 


Jim  *<*>?)  -  5(xo)  +  ffi  CxJ  JO  -  £j(x0jO  , 
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We  let  ^  be  the  unit  vector  of  the  direction  in  which  the  wire  points  outward  from  the  origin, 
and  we  let  /0  be  the  current  flow  out  of  the  origin.  Then  (2.36)  and  (2.37)  become 
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It  is  seen  that  L\(xq)  is  a  generalization  of  ATiOco)  as  defined  by  (2.31);  if  we  define 
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then  (2.39)  beconriu 
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which,  since  by  (2.43)  and  by  (2.44) 
fc*Po-0  , 
becomes 
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Reexpressing  (2.35)  with  an  eye  toward  symmetry,  we  have 
^^SCfo+yJ-ll^  +  jlSOcoiyJ  +  ^iCxo^  +  ^ijOfo^)  ,  (2.58) 

where  the  perturbation  terms  symmetric  with  respect  to  simultaneous  reversal  of  la  and  #  are 
Halo 


*«J>- 


s>  x  y  2(fr  x  xoXpo-y) 

Pi 


P$ 


(2.59) 


and  those  that  are  antisymmetric  are 
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Taking  derivatives,  we  And  that 
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whereas 
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It  is  apparent  that  the  field  derivatives  associated  with  are  artifact-contributed.  On  the 
other  hand,  ^42  is  not  implicated  as  artifact-producing  by  the  curl  test,  and  is  eliminated  on  the 
grounds  that,  like  Bs  and  unlike  B^\,  it  becomes  singular  everywhere  along  the  axis  of  the 
wire.  Hence,  in  (2.60)  we  have  an  expression  that  can  be  used  for  the  correction  of  artifact 
effects  in  the  elements  of  the  tensor  gradient  of  the  magnetic  field.  It  is  undesirable  to  make 
an  artifact  correction  to  the  magnetic  field  itself,  since  in  the  nonphysical  truncation  of  a 
current  flow  one  just  loses  the  field  contribution  from  the  truncated  section  of  the  current  flow 
without  also  gaining  a  nonphysical  field  component  in  its  place;  it  is  essentially  the  absence  of  a 
field  component  where  one  would  expect  it  on  physical  grounds  to  be  that  is  responsible  for  the 
nonphysical  behavior  of  the  field  derivatives  in  the  cose  of  the  semi-infinite  wire. 
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2.3  Streamline  tracing 

In  numerical  exploration  work  it  is  frequently  desirable  to  trace  streamlines  of  electri¬ 
cal  current  flow,  perhaps  to  check  for  peculiarities  in  the  flow  distribution  that  might  require 
special  programming  to  handle.  The  simplest  approach  to  streamline  tracing  is  to  pick  a  starting 
point,  find  the  current  density  vector  at  this  point,  extrapolate  in  the  direction  of  the  vector 
from  the  starting  point  over  some  convenient  interval  of  distance  to  a  new  point,  and,,  using 
this  new  point  as  a  starting  point,  repeat  the  process  for  as  tong  along  the  streamline  as  on 
wishes  to  trace;  however,  in  problems  where  there  is  some  region  in  which  the  direction  of  the 
current  density  vector  changes  rapidly  with  position  (a  characteristic  of  most  problems  of  practi¬ 
cal  interest),  this  method  can  lead  to  large  cumulative  errors  unless  the  step  size  is  kept  expen¬ 
sively  smalt  To  deal  with  this  difficulty,  a  more  complicated  method  of  the,  predictor-corrector 
type  was  developed,  and  was  found  to  be  satisfactory. 

We  assume  two  parallel  planes  separated  by  a  distance  d  and  a  streamline  passing 
through  them  both,  intersecting  the  first  plane  at  the  point  jco  and  then  the  second  plane  at  the 
point  Xl\  the  electrical  current  density  is  given  for  everywhere  between  the  planes.  The  basic 
idea  of  the  method  is  that  if  the  direction  of  the  current  density  vector  at  the  point  X£  along 
the  streamline  midway  between  the  two  planes  Is  determined  and  a  line  is  extended  from  Xo  in 
the  direction  of  »£  until  it  intersects  the  second  plane,  then  the  point  of  intersection  will  be 
very  close  to  JR,  provided  that,  the  curve  of  the  streamline  is  roughly  uniform  and  reasonably 
gentle.  To  approximate  the  position  of  this  midpoint,  we  first  determine  the  current  density 
vector  vo  at  the  point  xo,  then  draw  a  tine  from  xo  in  the  direction  of  vo  until  it  Intersects  the 
second  plane  at  the  point  xj; 


xi  —  d  +  xo  ,  (2.64) 

flv0 

where  Hist  unit  vector  perpendicular  to  the  planes  and  pointing  from  the  first  plane  toward 
the  second.  This  tine  will  miss  XI  on  the  convex  side  of  the  streamline.  We  then  find  the  point 
xj  midway  between  xo  andjq, 

xi-  xh  (xr-xj)  +  xi  (2.65) 
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and  find  v2  the  current  density  vector  at  x2;  this  vector  should  be  u  somewhat  better  approxima¬ 
tion  to  than  is  vjj.  We  then  draw  a  line  from  xq  in  the  direction  of  ^  until  it  intersects  the 
second  plane  at  the  point  xj; 
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Xj  —  d  ‘-*Lr  +  Xo  . 


(2.66) 


This  line  should  miss  jK  on  the  concave  side  of  the  streamline.  The  point  xj  midway  between 
xj  and  xo  is  given  by 


xj  -  VS  (xj-xj)  +  xj 
2  fl'va 

The  point  Xj  midway  between  xj  and  xj, 
jcj—  V6  Cxj-xj)  +  xj 
-  Vi  (xj+ xj)  , 
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should  then  be  fairly  dose  to  ,  and  the  current  density  vector  vj  at  xj  should  be  a  good 
approximation  to  t£.  Finally,  we  draw  a  line  from  xj  in  the  direction  of  vj  until  it  Intersects 
the  second  plane  at  the  point  xj; 
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The  point  xj  within  the  second  plane  should  then  lie  very  close  to  Xt,  Although  this  method 
requires  considerably  more  calculation  per  step  than  does  the  simple  extrapolation  from  xq 
along  the  direction  of  vj,  it  permits  one  to  safely  Increase  the  step  size  to  the  point  where  the 
total  amount  of  calculation  needed  to  trace  a  given  length  of  streamline  is  much  less  than  that 
required  by  the  simple  extrapolation  method. 

For  problems  in  which  a  streamline  may  double  back  upon  itself,  giving  rise  to  the 
possibility  of  denominators  going  to  zero  in  some  of  the  expressions  given  above,  this  method 
may  be  modified  to  sujt.  In  this  case,  we  picture  a  sphere  of  radius  d  centered  at  the  point  xj, 
and  postulate  a  streamline  passing  through  xj  and  emerging  through  the  surface  of  the  sphere 
at  the  point  XI,  The  formulas  corresponding  to  (2.64)  through  (2.69)  are 
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where  xj  b  the  method's  approximation  to  xo.  Note  that,  due  to  scaling  opportunities  not 
present  in  the  previous  case,  xl  and  x!J  in  (2.74)  may  be  replaced  with  x\  and  xl  respectively, 
thereby  eliminating  the  need  to  calculate  xj. 
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3.  Magnetic  anomalies  in  the  air  above  a  submerged  hemisphere 

In  order  to  gain  some  feel  for  the  orders  of  magnitude  of  the  magnetic  leld  values 
associated  with  a  localized  bottom  feature,  we  consider  the  case  of  a  hemisphere  with  its  base 
resting  on  an  otherwise  level  ocean  floor  and  with  its  radius  equal  to  half  of  the  ocean  depth. 
An  electrics]  current  density  field,  uniform  at  large  lateral  distances  from  the  hemisphere,  is 
induced  within  the  ocean  by  an  outside  source  of  EMF,  and  the  anomalous  magnetic  effects 
associated  with  the  bottom  feature  are  defined  to  be  those  arising  from  the  perturbation  on  the 
uniform  electrical  current  field  caused  by  the  feature. 

The  calculation  of  the  electrical  current  field  is  based  on  the  placement  of  current 
dipoles  in  such  a  way  as  to  satisfy  the  boundary  conditions  of  the  problem,  these  being  that 
there  should  be  no  current  flow  through  either  the  ocean  floor  or  the  ocean  surface.  The  basic 
formulae  for  the  estimation  of  the  appropriate  dipole  strengths  are  those  for  the  electrical 
current  density  distribution  about  and  within  a  sphere  of  radius  a  suspended  in  an  infinite 
ocean  within  which  the  electrical  current  density  field  approaches  Jok  at  large  distances  from 
the  sphere, 


7(rj)  -  Jo  I  k  -  corf  +  ‘/i  $  slnfl)  J  ,  (3.1) 

7(rjB)-0  ,  r<a  ,  (3.2) 

where  r  is  the  distance  from  the  center  of  the  sphere  and  the  direction  of  0-0  is  parallel,  to  A; 
expression  (3.1)  is  the  gradient  of  the  potential 

-  /0  J  f  +  ^  cosfl  *  (3.3) 

which  is  a  solution  to  Laplace’s  equation,  and  (3.2)  derives  from  a  trivial  solution  to  ^place's 
equation.  It  may  be  verified  that  on  the  spherical  surface  defined  by  r~a  the  curren  density 
field  component  J’t  normal  to  the  surface  is  equal  to  zero.  Expression  (3,1)  is  seen  n  be  the 
superposition  of  a  uniform  field  Jok  and  a  dipole  field  having  a  strength  of  and  a 

direction  antiparallel  to  the  uniform  component.  These  expressions  also  apply  to  the  case  of  a 
hemisphere  of  radius  a  resting  on  the  floor  of  an  ocean  of  arbitrarily  great  depth,  as  the 
cylindrical  symmetry  of  the  sphere  problem  about  the  x-axis  guarantees  that  no  streamline  will 
intersect  a  piano  containing  the  x-axis.  The  cose  of  a  hemisphere  resting  on  the  floor  of  an 
ocean  of  finite  depth  may  be  approximately  dealt  with  by  setting  up  an  array  of  many  dipoles  of 
equal  strength  and  a  common  orientation  antiparallel  to  the  uniform  field  Jok,  equally  spaced 
along  a  line  perpendicular  to  Jr  with  a  spacing  equal  to  twice  the  depth  of  the  ocean.  If  the 


ocean  floor  is  taken  to  correspond  to  the  plane  bisecting  this  line  at  the  central  dipole  of  this 
array,  and  the  ocean  surface  to  the  parallel  plane  midway  between  this  dipole  and  the  dipole 
immediately  above  it,  then  by  symmetry  there  is  no  current  flow  through  either  ocean  bottom 
or  ow*  urface;  however,  due  to  the  effects  of  the  dipoles  on  either  side  of  the  central  dipole, 
the  surface  about  this  dipole  on  which  7-f-O  is  no  longer  a  perfect  sphere,  but  an  oblate 
spheroid.  The  degree  of  this  hemispherical  distortion  was  examined  by  use  of  a  streamline  trac¬ 
ing  program  with  graphical  output,  and  it  was  found  for  an  array  of  five  dipoles  that  if  the  value 
of  a  is  no  more  than  half  of  the  ocean  depth  then  the  deviation  of  the  surface  about  the  central 
dipole  from  a  hemispherical  shape,  is  difficult  to  see  with  the  unaided,  eye,  implying  a  distortion 
of  no  more  than  about  5%  of  a , 

Granted  that  this  method  adequately  gives  the  electrical  current  density  distribution 
within  the  ocean,  the  next  step  in  the  magnetic  anomaly  calculations  is  to  use  expressions 
(2.10)  and  (2.11),  giving  the  magnetic  effects  due  to  the  current  density  distribution  within  a 
rectangular  box,  for  incremental  integration  over  the  current  density  field  to  estimate  the  total 
magnetic  field  and  its  spatial  derivatives  at  various  points  of  interest.  The  procedure  thut  was 
constructed  to  do  this  starts  with  . the  selection  of  a  volume  of  integration  having  a  rectangular 
horizontal  cross-section  and  reaching  vertically  from  the  ocean  floor  to  the  surface*  covering  a 
substantial  area  of  the  ocean  floor  about  the  bottom  feature;  this  volume  is  oriented  so  that  its 
vertical  faces  are  either  parallel  or  perpendicular  to  the  Si  direction.  The  volume  is  then  sec¬ 
tioned  by  a  set  of  equally  spaced  planes  perpendicular  to  the  Si  direction,  and  also  by  another 
set  of  equally  spaced  planes  (not  necessarily  the  same  spacing  as  the  first  set)  both  vertical  and 
perpendicular  to  the  first  set,  thus  breaking  the  horizontal  cross  section  of  the  volume  up  into  a 
regular  rectangular  grid.  Then,  one  goes  to  the  bottom  edge  of  the  upeurrent  face  of  the 
volume  of  integration,  and  draws  a  horizontal  line  along  the  face  slightly  above  the  edge,  mark¬ 
ing  the  points  along  this  line  that  Intersect  the  sectioning  planes  parallel  to  Jr;  these  points  serve 
as  starting  points  for  the  tracing  of  a  streamline  surface,  the  algorithm  for  which  is  given  in  the 
next  paragraph.  The  objective  of  this  part  of  the  procedure  is  to  trace  out  a  streamline  surface 
running  very  close  to  the  ocean  bottom,  and  use  this  surface  in  subsequent  calculations  as  the 
base  of  the  volume  of  integration;  this  is  substantially  less  cumbersome  than  numerically  speci¬ 
fying  the  shape  of  the  ocean  bottom,  and  the  error  introduced  into  the  calculations  by  this 
approximation  may  be  made  negligible  by  starting  the  streamline  surface  tracing  algorithm 
sufficiently  close  to  the  actual  bottom.  The  tracing  algorithm  yields  the  points  of  intersection 
between  the  streamline  surface  and  the  regular  array  of  vertical  lines  which  are  defined  by  the 
intersection  of  the  two  perpendicular  sets  of  sectioning  planes.  One  constructs  a  rectangular 
volume  about  each  one  of  these  vertical  tines  from  the  point  where  it  Intersects  the  ocean 


bottom  to  the  point  where  it  intersects  the  ocean  surface,  and  then  then  divides  this  into  small 
rectangular  boxes,  to  each  of  which  the  expressions  (2.10)  and  (2.11)  are  applied  and  the  mag¬ 
netic  field  contributions  therefrom  summed;  the  lateral  dimensions  of  the  rectangular  volume 
taken  to  be  about  each  line  are  determined  by  the  spacing  of  the  lines,  and  chosen  so  that  the 
vertical  faces  of  these  volumes  meet,  thereby  covering  the  whole  volume  of  integration.  The 
current  density  distribution  used  with  these  expressions  is  the  anomalous  current  density  distri¬ 
bution,  which  is  the  total  field  with  the  uniform  component  J<#  subtracted  off;  the  current  den¬ 
sity  field  between  the  ocean  surface  and  the  streamline  surface  is  then  only  the  superposition  of 
the  individual  dipole  fields,  and  the  field  below  the  streamline  surface  and  above  the  ocean  floor 
is  just  -Jro&»  with  the  field  everywhere  else  being  equal  to  zero. 

The  algorithm  used  for  tracing  a  streamline  surface  is  partly  based  on  the  first  stream¬ 
line  tracing  algorithm  discussed  in  section  2.2.  We  suppose  that  the  sectioning  planes  perpendic¬ 
ular  to  £  are  separated  front  each  ether  by  a  distance  d,  and  that  on  one  of  them  we  have  a  set 
of  points  along  the  intersection  curve  of  the  plane  with  the  streamline  surface,  these  points 
being  the  points  of  intersection  between  this  curve  and  the  vertical  tinea  defined  by  the  inter¬ 
section  of  the  plane  and  the  perpendicular  set  of  sectioning  planes.  Using  a  single  step  of  the 
streamline  tracing  algorithm  for  each  one  of  these  points  taken  as  a  starting  point,  we  get 
another  set  of  points  on  the  next  plane  over  in  the  *  direction,  to  which  we  can  fit  another, 
slightly  different,  curve;  in  general,  these  points  will  not  lie  on  the  vertical  lines  of  intersection 
in  this  plane.  Finally,  we  go  in  turn  to  each  of  the  vertical  lines  of  intersection  of  the  second 
plane,  and  fit  a  cubic  polynomial  curve  through  the  four  points  nearest  it  (those  usually  being 
the  two  nearest  points  on  either  side  of  the  line),  and  and  record  the  point  where  the  fitted 
curve  crosses  the  line.  This  last  set  of  points  may  be  used  to  start  another  step  of  this  algo¬ 
rithm,  which  may  be  used  recursively  to  trace  the  points  of  the  entire  streamline  surface.  As  a 
final  check  of  accuracy,  one  may  note  whether  or  not  the  streamline  surface  tends  toward  a  flat 
sheet  slightly  above  the  level  ocean  floor  as  it  is  extended  past  the  localized  bottom  feature. 
Note  that  it  is  unnecessary  to  trace  out  and  store  the  points  for  the  entire  streamline  surface 
before  proceeding  with  the  integration;  it  is  possible  to  perform  the  integration  for  a  given 
plane  after  computing  the  streamline  surface  points  of  that  plane  from  those  of  the  previous 
plane,  add  the  results  of  the  integration  to  an  accumulated  sum,  throw  away  the  points  of  the 
previous  plane,  and  then  go  on  to  the  next  plane,  and  so  on,  with  considerable  savings  in  pro¬ 
gram  array  storage  space. 

In  dealing  with  the  electrical  current  density  distribution  truncation  error  problem  dis¬ 
cussed  in  section  2.3,  a  method  that  has  proved  reasonably  effective  in  practice  is  to  examine, 
in  the  course  of  integration,  each  of  the  rectangular  boxes  that  constitute  the  busic  unit  of 
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integration  to  see  ir  it  has  one  or  more  faces  on  the  outside  of  the  volume  of  integration;  if  an 
outside  face  is  detected,  then  the  current  density  vector  at  its  center  is  multiplied  by  the  area  of 
the  face  in  order  to  estimate  the  total  current  flow  through  the  face,  and  then  equation  (2.60) 
from  section  2.3  is  used  with  this  information  to  calculate  corrections  to  the  elements  of  the 
tensor  gradient  of  the  magnetic  held.  Failure  to  make  such  corrections  will  result  in  nonzero 
values  for  the  curl  of  the  magnetic  field  in  places,  such  as  empty  air,  where  there  Is  no  electrical 
current;  practical  experience  has  shown  that  these  tensor  element  corrections  can  be  comparable 
in  magnitude  to  the  elements  of  the  corrected  tensor.  A  further  problem  exists  in  that  the 
current  density  held  on  the  faces  of  a  box  is  estimated  by  the  imperfect  method  of  extrapola¬ 
tion  from  its  value  and  tensor  gradient  at  the  center  of  the  box,  and  that  therefore  one  can 
expect  slight  discontinuities  in  the  current  density  field  as  one  passes  from  one  box  to  the  next;  ' 
properly,  the  effects  of  these  discontinuities  should  be  corrected  for  with  further  applications  of 
(2.60).  Also,  there  will  be  minor  errors  from  the  use  of  (2.60)  on  the  outside  box  faces  with 
what  amounts  to  face-averaged  current  densities,  rather  than  the  more  correct  procedure  of 
Integrating  (2.60)  over  the  surfaces  of  the  faces  in  question.  However,  both  of  these  additional 
sources  of  error  can  be  reduced  by  cutting  down  on  the  size  of  the  basic  unit  of  integration, 
and,  with  the  primary  source  of  error  allowed  for  as  described,  practical  calculations  imply  that 
accuracy  to  two  significant  places  seems  feasible. 

The  actual  calculations  were  done  for  a  roughly  hemispherical  shape  of  radius  4.0 
meters  at  its  base  on  the  bottom  of  an  ocean  4.3  meters  deep  (see  figure  2.1),  assuming  a  uni¬ 
form  electrical  current  density  field  component  of  1.0  ampere/ meter2  in  the  %  direction;  scaling 
relationships  may  be  used  with  the  fields  and  their  tensor  gradients  calculated  for  this  model  to 
estimate  the  corresponding  values  for  a  somewhat  larger  system  with  a  different  uniform 
current  density  field  component.  The  flattening  shown  in  the  figure  near  the  top  of  the  projec¬ 
tion  from  the  seafloor  is  the  result  of  the  use  of  only  five  dipoles  to  simulate  the  electrical 
current  flow  about  a  hemisphere  in  the  close  vicinity  of  the  plane  of  symmetry  which  in  this 
representation  stands  for  the  ocean  surface.  Had  a  much  larger  number  of  dipoles  been  used  in 
the  simulation,  this  flattening  effect  could  have  been  greatly  reduced,  yielding  the  flow  pattern 
about  a  nearly  perfect  hemispherical  surface;  however,  such  a  flattening  at  the  top  is  in  fact  a 
characteristic  of  some  oceanic  features  of  interest,  such  ns  seamounts,  and  accordingly  is  more 
realistic  for  features  of  this  sort.  The  basic  unit  of  integration  was  a  cubical  box  10  centimeters 
on  a  side,  and  the  grid  spacing  of  the  vertical  lines  of  intersection  of  the  sectioning  planes  was 
10  centimeters  by  10  centimeters.  At  each  one  of  these  vertical  lines,  cubes  were  stacked  from 
the  ocean  floor  upward  until  a  cube  extended  beyond  the  streamline  surfuce,  at  which  point  this 
cube  was  truncated  level  with  the  streamline  surfuce;  expressions  (2,10)  and  (2,11)  were  then 


FIGURE  2,1 


Streamline  flow  over  the  roughly  hemispherical  bottom  feature  in  the  vertical  plane  parallel  to 
the  flow  at  a  distance  from  the  feature  and  passing  through  the  centerline  of  the  feature;  to 
point  up  the  extent  of  the  deformation  from  the  flow  pattern  over  a  perfect  hemisphere,  a  hem¬ 
isphere  of  the  appropriate  radius  is  superposed  over  the  streamline  pattern. 
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used,  assuming  a  uniform  current  density  field  of  -1.0  ampere! meter1  x  within  all  of  the 
integration  units.  Then,  cubes  were  stacked  from  the  top  of  the  truncated  cube  until  one 
extended  above  the  ocean  surface,  at  which  point  this  cube  was  truncated  levei  with  the  sur¬ 
face;  the  current  density  fields  and  their  tensor  gradients  used  in  expressions  (2.10)  and  (2  11) 
are  generated  by  the  five  dipoles  making  up  the  vertical  dipole  array  for  this  problem.  The  total 
volume  of  integration  was  20  meters  by  20  meters  laterally  by  4.5  meters  vertically,  with  the 
hemisphere  located  at  its  bottom  center.  Tracing  of  the  streamline  surface  was  started  15  cen¬ 
timeters  over  the  ocean  bottom  at  the  upcurrent  edge  of  the  volume  of  integration,  and  its 
extrapolation  at. the  downcurrent  edge  of  this  volume  differed  no  more  than  5  centimeters  from 
this  height  at  any  point  along  the  upcurrent  side,  verifying  that  the  extrapolation  hta  been 
properly  done.  Calculations  were  done  for  various  points  of  interest  in  the  air  over  the  ocean 
surface,  and  the  results  are  presented  in  table  2.1.  The  x-axis  of  the  coordinate  system  points  in 
the  direction  of  the  uniform  component  of  the  electrical  current  flow,  and  the  z-axis  points 
upward;  the  origin  of  the  coordinate  system  is  on  the  ocean  bottom,  at  the  center  of  the  hemi¬ 
sphere. 


We  make  this  work  useful  in  practice  by  developing  scaling  relations  for  the  magnetic 
field  and  its  tensor  gradients  in  going  from  one  system  with  one  length  scale  to  another  system 
with  a  different  length  scale  but  with  the  same  ratios  between  its  various  length  parameters  as 
the  first  system.  Referring  to  equation  (2.1),  we  see  that  the  magnetic  field  is  linear  with 
respect  to  the  electrical  current  density.  Further  work  with  this  equation,  experimenting  with 
changes  of  length  variables  by  multiplicative  factors,  indicate  that  if  we  have  two  systems, 
identical  except  as  to  scale,  and  By,  Ju  and  d\  are  magnetic  field,  electrical  current  density,  and 
length  parameters  of  the  first  system  whereas  B J2,  and  a2  are  the  corresponding  parameters 
of  the  second  system,  then  we  have 


By  ./j  d  j 

2?i  J\  d\ 


(3.4) 


the  same  process  applied  to  spatial  derivatives  of  equation  (2.1)  shows  that  if  By  arid  B{  ar: 
corresponding  magnetic  field  gradients  for  the  two  systems,  then 


B£  h 
By  “  J\ 


(3.5) 


which  is  to  be  expected,  as  the  larger  system  will  show  a  smaller  percentage  change  in  a  field 
value  than  the  smaller  system  for  a  given  change  in  position. 

As  an  example,  suppose  that  we  wish  to  find  the  magnetic  field  and  its  tensor  gradients 
at  a  point  1.0  kilometer  in  the  air  over  the  surface  of  an  ocean  4,5  kilometers  deep,  with  a 
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y  — 4.0  meters  x  — 0.0  meters  x<™  4.0  meters 


B, 

0.000 

-0.111 

By 

0.132 

0.137 

B 

-0.338 

-0.194 

B» 

-0.048 

0.002 

By 

0.000 

•0.001 

Be 

0.000 

0.048 

Byy 

-0.084 

-0.055 

By,  • 

0.016 

-0.007 

Bn 

0.132 

0.053 

y  -  0.0  meters 

X'  —  0,0  meters 

x  —  4.0  meters 

Be 

0.000 

0.000 

B , 

0.473 

0.325 

B 

0.000 

0.000 

Bex 

0.000 

0.000 

Byy 

0.000 

-0.050 

Be 

0.000 

0.000 

Byy 

0.000 

0.000 

By, 

-0.181 

-6.088 

Bn 

0.000 

0.000 

TABLE  2.1 

Field  and  field  gradient  values  at  selected  points  in  a  horizontal  plane  1.0  meter  above  the  sur¬ 
face  of  an  ocean  4.S  meters  deep  with  a  hemispherical  bulge  of  radius  4.0  meters  resting  on  the 
bottom,  as  described  in  the  text  on  page  36,  The  field  values  are  in  unitr  J  10 ~6 tesla,  and  the 
field  gradient  values  are  in  units  of  1  Or" tesla /meter.  The  origin  of  the  plane  is  on  the  axis  of 
the  bulge.  In  the  absence  of  the  bulge  the  field  would  be  a  uniform  B--2.82X  lGr6j>  tesla . 


shape  such  as  illustrated  in  figure  2,1,  having  a  radius  at  its  hase  of  4.0  kilometer:.  resting  on 
the  ocean  floor  directly  below  the  point;  we  assume  a  uniform  electrical  conductivity  ui'  3.3  S/m 
for  the  ocean,  and  suppose  that  sources  of  EMF  outside  of  the  ocean  induce  an  essentially  con¬ 
stant  electric  field  within  it  that  is  a  uniform  1,0  n  V/tri  in  the  x  direction  at  great  distances 
from  the  distorting  effects  of  the  hemisphere.  The  system  described  is  that  of  table  2.1  scaled 
up  in  size  by  a  factor  of  1000,  With  an  electrical  current  density  at  large  distances  from  the 
hemisphere  of  3.3x10"*  ampere /meter1  in  the  k  direction,  as  opposed  to  1.0  ampere/ meter1  for 
the  smaller  system.  From  equation  (3.4)  we  then  have  for  the  field  values  that 
3.3xl0"35i,  and  from  equation  f3.5)  we  have  for  the  field  value  gradients  that 
Bi'-  3.3x  10"*  B\.  Referring  to  table  2.1  to  get  the  field  values  directly  above  the  center  of 
the  bottom  feature,  we  therefore  have  for  By,  the  only  nonzero  fiold  component,  a  value  of 
1.6xl(T9  tesla,  and  for  By,  and  f?v,  the  only  nonzero  tensor  gradient  elements,  a  common 
value  of  — 6.0x  10"13  tesh  per  meter.  Using  this  procedure,  we  scale  up  the  contents  of  table 
2.1,  giving  the  results  in  table  2.2. 

As  there  are  np  electrical  currents  in  the  air  above  the  ocean  surface,  the  magnetic 
field  in  this  region  can  be  expressed  as  the  gradient  of  a  scalar  potential  which  is  a  solution  to 
Laplace’s  equation.  The  principle  restriction  on  the  practical  use  of  the  method  of  magnetic 
field  computation  discussed  in  thtTtext  above,  that  the  points  at  which  the  field  is  computed 
must  be  a  distance  above  the  ocean  surface  that  is  on  the  order  of  many  times  the  largest  linear 
dimension  of  the  basic  unit  volume  of  integration,  can  in  principle  be  circumvented  by  fitting  a 
general  series  expansion  solution  to  Laplace’s  equation  to  the  values  computed  by  this  method, 
and  then  using  the  resulting  potential  expression  to  compute  the  magnetic  field  and  its  tensor 
derivatives  at  points  elsewhere;  probably  the  best  such  solution  to  Laplace  s  equation  for  this 
purpose  is  an  expansion  in  the  spherical  harmonics  about  the  center  of  the  base  of  the  hemis¬ 
pherical  shape,  using  only  those  terms  which  vanish  at  arbitrarily  large  radial  distances  from  this 
point.  However,  one  should  bear  in  mind  that  downward  continuations  of  this  sort  tend  lo  be 
numerically  unstable,  in  that  small  errors  in  the  determination  of  the  field  values  at  the  1.0 
kilometer  altitude  above  the  ocean  surface  can  easily  be  multiplied  to  large  errors  in  the  values 
computed  at  the  surface  from  the  potential  expression,  and  that  the  field  gradients  computed  by 
the  method  of  this  section  are  mildly  prone  lo  distortion  by  the  effects  of  those  mathematical 
artifacts  which  cannot  be  completely  compensated  (principally  those  introduced  at  the  faces  of 
the  basic  rectangular  volumes  of  integration  making  up  the  simulation  volume,  and  in  particular 
those  faces  which  are  on  the  outer  surface  of  this  volume);  it  could  easily  happen  in  many  cases 
that  the  continuation  error  exceeds  the  magnitude  of  the  value  being  computed,  and  it  is 
recommended  thut  in  fitting  the  potential  expression  to  the  values  numerically  determined  in 
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4.0  km 

x- 0.0  km 

x«4.0fcm 

*x 

0.00 

•0.37 

By 

0.44 

0.45 

Bt 

-U 

-0.64  . 

B» 

-0.16 

0.005 

B*> 

0.00 

-0.003 

Bn 

0.00 

0.16 

Byy 

By* 

-0.28 

•0.18 

0.05 

*0.02 

0.44 

0.1S 

y~0.0km 

X*  0.0  km 

x-4.0  km 

B* 

0.00 

0.00 

By 

1.6 

0.11 

B> 

0.00 

0.00 

Bxx 

0.00 

0.00 

By 

0.00 

-0,15 

B„ 

0.00 

0.00 

Byy 

0.00 

0.00 

By. r 

-0.60 

-0.29 

Bn 

0.00 

0.00 

TABLE  2.2 

Field  and  field  gradient  values  at  selected  points  in  a  horizontal  plane  1.0  kilometer  above  the 
surface  of  an  ocean  4,5  kilometers  deep  with  a  hemispherical  bulge  of  radius  4.0  kilometers 
resting  on  the  bottom,  as  described  in  the  text  on  pages  38  and  40.  The  held  values  arc  in  units 
of  10"9r«to,  and  the  field  gradient  values  are  in  units  of  10“ 12  tesla/ meter .  The  origin  of  the 
plane  is  on  the  axis  of  the  bulge.  In  the  absence  of  the  bulge  the  field  would  be  a  uniform 
B--9.3X  1< tesla. 


moderate  to  large  distances  above  the  ocean  surface  one  relies  heavily  on  the  magnetic  field 
vector  elements,  which  should  be  somewhat  less  contaminated  by  mathematical  artifacts  than 
the  tensor  gradient  elements. 


4.  Solution  for  the  electrical  current  density  distribution  in  an  ocean  with  arbitrary  bottom 
features 

Given  the  methods  of  section  2,  the  magnetic  anomalies  resulting  from  a  given  bottom 
feature  may  be  computed  provided  that  the  electrical  current  disturbances  associated  with  that 
feature  are  known.  In  this  section  an  algorithm  is  given  for  finding  the  electrical  current  distri¬ 
bution  about  a  localized  bottom  feature  of  arbitrary  shape. 

The  construction  process  starts  with  a  level-surfaced  ocean  of  indefinite  depth,  within 
which  them  is  a  horizontal  electric  field  presumed  to  be  induced  by  sources  in  the  upper  atmo¬ 
sphere.  This  primary  electric  field  £o  is  taken  to  be  the  gradient  of  the  potential  </>0.  «nd  is 
assumed  to  induce,  an  electrical  current  distribution  Jo  according  to  the  relation 

JJ-o-V^o  .  (4.1) 

where  cr  is  the  electrical  conductivity  distribution  for  the  ocean. 

A  pillbox-shaped  closed  surface  is  then  postulated,  whose  top  is  the  level  ocean  sur¬ 
face  and  whose  bottom  has  the  shape  of  a  hypothetical  ocean  floor  profile;  this  profile  is  taken 
to  be  level  and  parallel  to  Jj>  except  for  a  region  close  tp  the  axis  of  the  pillbox,  where  the 
profile  is  allowed  to  take  any  shape  that  does  not  intersect  the  top  of  the  pillbox.  We  now 
assume  a  potential  <£1  and  an  electrical  current  density  distribution  7  such  that 

7-aV^|)  ,  (4.2) 

and  require  as  boundary  conditions  on  7  that  its  component  nbrmal  to  the  pillbox  surface 
should  be  zero  at  all  points  of  the  top  and  bottom  of  the  box  (indicating  the  absence  of  current 
flow  through  either  boundary),  and  that  the  pillbox  sides  should  be  far  enough  removed  from 
the  central  disturbance  that  7  may  safely  be  assumed  to  be  equal  to  T0. 

The  boundary  condition  on  7  at  the  top  of  the  box  is  already  met  by  T0,  so  we  have 
for  0]  there  that 

5-V0i  -  0  ,  (4.3) 

and  at  the  bottom  of  the  box  we  immediately  have  from  (4.2)  that 

A-V0o  .  (4.4) 

where  h  is  the  inward-pointing  unit  vector  normal  to  the  box  surface;  since  by  hypothesis  Jo  is 
parallel  to  the  bottom  except  near  the  axis  of  the  box,  this  reduces  to  (4.3)  over  most  of  the 
box’s  bottom.  On  the  box  sides,  the  condition  that  7- Jo  implies  by  (4.1)  and  (4.2)  that 


The  normal  derivative  of  <f>\  is  therefore  specified  over  the  entire  closed  surface  of  the  pillbox, 
thereby  constituting  Neuman  boundary  conditions,  which  are  sufficient  to  define  a  solution  for 
unique  to  within  an  additive  constant.  We  thus  have  7  everywhere  within  the  box;  as  7 
meets  the  boundary  conditions  appropriate  to  an  electrical  current  density  distribution  within  an 
ocean  having  the  specified  bottom  topography,  the  problem  is  solved.  For  methods  to  compute 
given  Neuman  boundary  conditions,  see  Forsythe  and  Wasow,  1960. 


S.  A  quick,  inexact  method  for  approximation  of  the  electrical  current  density  distribution 
due  to  a  localized  bottom  feature 

For  situations  where  the  aspect  of  the  bottom  feature  is  small,  a  quick  if  inexact 
method  for  estimating  the  associated  electrical  current  flow  disturbance  is  outlined.  As  in  the 
previous  section,  one  assumes  a  horizontal  electric  field  such  as  would  exist  in  the  absence 
of  an  ocean  bottom,  constructs  a  hypothetical  surface  within  the  ocean  having  the  shape  of  the 
actual  bottom  topography,  then  calculates  at  each  point  on  this  sheet  the  component  of  Eq 
parallel  to  the  surface  normal.  A  surface  charge  distribution  over  the  bottom  sheet  is  then 
estimated  that  will  give  rise  to  an  electric  field  whose  normal  component  at  each  point  on  the 
sheet  exactly  cancels  the  corresponding  normal  component  of  fij.  The  supplementary  condition 
that  there  be  no  electrical  Current  flow  through  the  ocean  surface  (which  is  assumed  for  con¬ 
venience  to  be  level)  may  be  met  by  constructing  a  system  of  image  reflections  of  the  estimated 
surface  charge  distribution  through  the  planes  of  the  ocean  surface  and  the  level  ocean  bottom 
from  which  the  bottom  features  are  taken  to  protrude.  This  method  is  suggested  by  the  physi¬ 
cal  observation  that  an  electrical  current  flow  over  a  bottom  irregularity  will  in  fact  deposit 
charges  on  the  irregularity’s  surface  in  a  fashion  tending  to  divert  the  flow  away  from  the  irre¬ 
gularity. 

.  There  appears  to  be  no  rapid,  straightforward  algorithm  for  deriving  a  surface  charge 
distribution  that  will  exactly  solve  thic  above  problem.  However,  an  approximate  method  exists 
which  should  serve  adequately  in  many  cases  of  practical  interest.  Suppose  that  we  have  an 
infinite  plane  surface  with  a  uniform  surface  charge  density  pu  ,  which  is  embedded  in  an 
infinite  uniform  dielectric  medium  with  a  dielectric  constant  then  by  Gauss’  law  we  have 
that  the  electric  field  in  the  medium  is  given  by  the  relation 

(5.1) 

where  ft  is  the  normal  unit  vector  outward  from  the  plane.  If  the  plane  suffers  a  local  deforma¬ 
tion  of  small  aspect,  then  this  should  still  be  a  reasonably  good  approximation  for  the  electric 
field  at  distances  from  the  surface  small  compared  to  the  lateral  extent  of  the  deformation. 

Similarly,  if  the  surface  charge  is  allowed  to  vary  gradually  from  one  point  to  another 
on  the  surface*  then  for  our  purposes  the  approximation  should  remain  serviceable  for  distunccs 
from  the  surface  small  with  respect  to  the  scale  of  variation;  such  a  charge  variation  would  pro¬ 
duce  a  horizontal  electric  field  component  close  to  the  surface,  but  the  local  field  indicated  by 
(5.1)  should  nevertheless  be  a  fairly  good  estimate  of  the  normal  component  of  the  field  at 
small  distances  from  the  surface,  The  seawater  medium  with  which  we  are  in  practice 
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concerned  is  not  a  dielectric  material,  but,  as  it  is  only  weakly  conducting  and  in  principle  the 
distribution  of  free  charges  in  the  problem  is  maintained  in  a  steady  state  by  outside  sources  of 
EMF  (upper  atmospheric  sources,  etc.),  it  may  be  considered  to  be  a  dielectric  material  for  the 
purposes  of  these  calculations.  The  effective  dielectric  constant  of  this  medium  is  difficult  to 
estimate,  and  it  is  fortunate  that  its  value  has  no  effect  On  the  end  results  of  the  calculations.  In 
calculations  deriving  the  surface  charge  p  from  the  held  Eq  within  a  medium  of  dielectric  con¬ 
stant  «  it  is  the  ratio  <p/«)  that  is  solved  for,  and  it  is  this  ratio  that  is  used  in  turn  to  compute 
the  correction  held  to  provided  that  the  surface  charge  is  used  only  as  an  intermediate  cal¬ 
culations!  device,  rather  than  being  given  any  physical  significance,  any  value  of  c  will  do  pro¬ 
vided  that  it  is  consistently  used.  A  value  of  90  *0>  roughly  the  value  for  distilled  water  ai  0°C, 
would  be  a  good  choice  if  one  wishes  p  to  have  a  physically  reasonable  range  of  values. 

The  recommended  approximation  is,  then,  given  by 
p--2e,£'0A  .  (5.2) 


where  h  is  the  normal  unit  vector  directed  outward  from  the  charged  surface  at  a  given  point, 
r0  is  the  uncorrected  held  vector  at  that  point,  and  «,  2s  the  chosen  value  for  the  effective 
dielectric  constant  of  the  seawater.  For  calculation  of  the  correction  held  E\  one  might  sweep 
the  surface  charge  distribution  into  a  number  of  small  piles  and  use  the  formula 


£i0 ?) 


h 


<ii 


(5.3) 


where  N  is  the  number  of  plies,  qi  is  the  accumulated  charge  in  each  pile,  and  iq  is  the  position 
of  each  pile.  If  the  surface  charge  distribution  derived  from  (5.2)  is  divided  sufficiently  finely, 
the  resultant  error  in  the  held  calculations  will  be  smaller  than  the  error  of  approximation  in 
using  (5.2)  to  get  the  surface  charge  distribution.  In  any  case,  it  is  strongjy  advised  to  use  a 
streamline  tracing  program  such  as  is  described  in  section  2  to  verify  that  the  above  procedure 
has  given  an  answer  acceptable  for  one’s  purposes,  In  the  event  that  the  approximate  solution 
falls  short  of  adequacy,  it  may  still  usefully  serve  as  a  starting  point  for  a  charge-redistributing 
variation  program  leading  to  a  more  exact  solution.  The  electric  current  density  distribution  7 
within  the  ocean  is  given  by 


7  ~<r  {£$¥£{)  , 


(5.4) 


where  <r  is  the  electrical  conductivity  distribution  within  the  ocean, 
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6.  Shallow  features  on  the  ocean  floor 


6.1  Basic  mathematics 

The  shape  of  any  reasonably  smooth  bottom  feature  can  be  described  by  the  superposi¬ 
tion  of  terms  of  the  form 

{( x%y)  -  a  expiiicH*)  ,  (6.1) 

£#  —  k  ,  k>0  ,  (6.2) 

where  £  is  altitude  relative  to  the  mean  depth  of  the  ocean  floor,  and  the  magnetic  effects  of  a 
single  such  tesrm  may  be  computed  by  means  of  a  variation  of  the  mathematics  developed  In 
section  3  of  chapter  1  for  finding  the  magnetic  effects  of  a  shallow  wave  on  the  ocean  surface. 
Provided  that  the  relation  between  bottom  profile  and  associated  magnetic  effect  is  linear,  the 
magnetic  field  due  to  a  bottom  profile  given  by  the  superposition  of  a  number  of  terms  of  the 
form  (6.1)  is  just  the  superposition  of  the  fields  calculated  for  the  individual  terms.  It  will  be 
seen  that,  for  this  linearity  to  be  safely  assumed,  the  amplitude  of  each  term  must  be  somewhat 
smaller  than  its  wavelength  and  somewhat  smaller  then  the  depth  of  the  ocean,  and  the  electri¬ 
cal  conductivity  of  the  seawater  at  the  ocean  bottom  must  be  essentially  constant  over  the 
term’s  amplitude  range. 

The  ocean  bottom  feature'  problem  differs  from  the  surface  wave  problem  discussed  in 
section  3  of  chapter  1  in  that  the  scale  of  the  induced  electrical  current  density  distribution  per¬ 
turbations  could  reasonably  he  a  significant  fraction  of  the  depth  of  the  ocean,  and  the  electrical 
conductivity  of  the  ocean  cannot  automatically  be  considered  uniform  throughout  the  ocean. 
Keeping  this  in  mind  In  the  adaptation  of  the  surface  wove  mathematics  to  the  present  prob¬ 
lem,  we  write  for  the  electrical  current  density  distribution  7 


f  —  cr  H!  (6.3) 

2?  «  V*  ;  (6.4) 

the  divergence  condition 

V-fl-Q  (6.5) 

then  gives  us  from  (6.4) 

Vfy  »  0  .  to  6) 


We  assume  an  electric  field  £q,  induced  by  sources  of  EMF  outside  of  the  ocean,  that  wo  ,.:d  be 


the  field  within  the  ocean  in  the  absence  of  bottom  features;  it  is  taken  to  be  horizontal  and 
uniform  within  the  ocean,  with  the  form 


Fq-Eo*  .  (6.7) 

We  then  have  for  the  general  solution  to  (6.6) 

(ft  —  £o*  +  A  exp  (Iki'll)  +  B  exp  (6.8) 

lc\- k  (cos9,sin0-i)  (6.9) 

£2- k  (cosBpinBJ)  ,  (6.10) 

and  hence 

E  -*  Erf  4-  £1  /  A  exp(/£l‘3?)  +  i  B  exp  •  (6.11) 


The  term  in  A  represents  a  field  component  that  grows  exponentially  upward,  and  the  term  in  B 
represents  1  component  that  damps  exponentially  upward;  if  the  ocean  were  assumed  to  be 
infinitely  deep,  then  the  term  in  A  would  vanish. 

Let  us  define  the  vertical  coordinate  so  that  the  mean  level  of  the  ocean  floor  is  at  z-0 
and  the  level  of  the  ocean  surface  is  at  z-d.  The  boundary  condition  at  the  surface,  that  there 
be  no  electrical  current  flow  through  the  level  surface,  is 

7-2—0  at  z—d  ,  (6.12) 

which  implies  through  (6.3)  that 

TS-2-0  at  z-d  ;  (6.13) 

from  this  conditio*  and  (6.11)  we  haVe  that 

A  -  JJ  exp  (-2kd)  .  (6.14) 


The  corresponding  boundary  condition  at  the  bottom  is 

7-7-0  ot  z-{Cx'iy)  , 


7- 


where  7  is  a  vector  normal  to  the  ocean  tioor  at 
from  (6.3)  we  have 


2“  i(x\y)  as  defined  by  (6.1). 


(6.15) 

(6.16) 
Once  again. 


1-7-  0  at  z-i  ; 
using  (6.11)  for  5 ,  assuming  that 
|to  |  «  1  , 
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(6.17) 


(6.18) 


and  discarding  terms  quadratic  or  higher  in'(ka),  we  get  from  (6.17) 
B  —  —  i  a  Eqx  coS0  , 

where 


(6.19) 


1 


*"  l- exp  (-2kd)  : 

We  then  have  from  (6.11)  for  the  electric  Held  within  the  ocean  subject  to  (6.18) 


(6.20) 


2?  -  £o  {•*  +  a  *t X  co.s0  txp (-2 kd)  exp (lk\‘T)  +  a  k\x  cost) exp (ik^-T)) 


(6.21) 


In  the  course  of  this  derivation  one  finds  that,  in  order  for  the  discarding  of  the  terms  higher 
than  linear  in  (ka)  to  be  valid,  one  must  also  assume  that  c!  >>  |£|.  In  particular,  If  a  >  xhd 
then  the  result  (6.19)  becomes  grossly  invalid. 

From  equations  (2.1)  and  (6.3)  we  have  for  the  magnetic  field  50c o)  at  the  position  x~o 


5(xo)  *■  '~rr  J  dx  J  <fy  f  it  <r(r) 

*WT  rM  -oa  ( 


5ce)  x  Oc^e) 


l*V*l3 


(6.22) 


where  the  electrical  conductivity  of  the  seawater  is  assumed  not  to  vary  laterally.  For  purposes 
of  calculation  we  find  it  expedient  to  separate  (6.22)  into  two  parts; 


5(*o)  -  5JOco)  +  5p(*5) 

*•  to  r  K  ,>*<*)*  OTr*) 

*  "  A  *  X  *  {  * 

Mo  ,,*«>*  (xtr%) 

So">- 47  /  *  /  *  4  *  ^  — UV5F 


(6.23) 

(6.24) 


(6.25) 


If  we  expand  the  integrand  of  (6,24) in  powers  of  z  about  z-0,  we  find  that  only  the  term  con¬ 
stant  with  respect  to  z  contributes' to  5J(x^)  less  than  quadratically  in  powers  of  (ka).  Accord¬ 
ingly,  retaining  only  this  term  and  letting  <r/j  be  the  value  of  ir(z)  at  z-0,  we  have  for  (6.24) 
to  first  order 


*♦/-*  HorgEo  f  j  f  ,  UVo ,\y-MI  \ 

S  Uo)  "  4 rr  £  dx  X  ^  i  (x--x0)l+ (y-.yo)1 X+*W  1  ^ 


(0,z0,(v-->'o)l 


(6,26) 
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Note  that  retaining  only  terms  up  to  first  order  in  (ka)  requires  us  to  disregard  the  effects  of 
variation  of  o-(z)  over  the  amplitude  range  of  the  bottom  feature.  Making  then  the  substitu¬ 
tions 


u  -  (x-xq)  cos®  +  (y-yo)  sin® 
v  -  (y-yo)  cose  -  (x-Xo)  sin® 


Uo-  x0  cos®  +  y0  sin®  , 


(6.27) 

(6.28) 
(6.29) 


we  get  from  (6.26) 


Moo’a^'ofl  V  r  .  r  .  tu,z0(u  sine  -t-  v  costrr  v 

"s<* o)  -  — 4^ — ■  expdkuo)  J.  du  dv  — - exp(/fcw)  ’  (6’30) 


(0,zo,u  sin®  +  v  cos®) 


from  which  integration  over  v  gives,  by  use  of  symmetry  and  equation  (3.30)  of  chapter  1, 


/r* WnEffl  /.,  s  r  ,  tu^o.M  sintr;  , 

Bs(xq)  -  — ^ —  exp(lku0)  J  du  — ui^.^ —  exp  Oku)  . 


(O^o sin®) 


(6.31) 


This  integral  yields  easily  to  contour  integration,  giving  for  z0>0  (the  only  region  of  practical 
interest) 


l5£(xo)  -  'huipBEva  expQkua)  exp(-kzo)  (0,1,/  sin®)  . 


(6.32) 


We  find  it  convenient  also  to  further  separate  By  into  the  parts 
2Tk(Xo)  Byu  (Xq),+  Byp\{Xo)  +  Byp jOcq)  , 


(6.33) 


where  Byu  is  concerned  with  the  uniform  component  of  (6.21),  ByP \  with  the  perturbation 
term  in  and  ByPJ  with  the  perturbation  term  in  Byu  has  already  been  evaluated  in  sec¬ 
tion  2  of  chapter  1.  For  B~yP\  we  have 


Jfc,«>  -  (-2W)  J  d*  f  dyt  <4 

W  _W  ^0 


x  <r(z)  expilkfi'X)  expikz) 


(cos6,$ln0—l)x  (xp-ff) 
ix>*|3 


(6.34) 


Using  the  substitutions  (6,27)  through  (6,29)  together  with  the  definitions 
0  “  k  cos fl  +  p  sin® 

D  -  P  cos®  -  x  sin®  , 


(6.35) 

(6.36) 


r,rf.rT 
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this  becomes 


Byr i(*o)“  ~M ^^'aL  X  C°—~  exp(-2kd)  exp  (iku0)  f  du  J  dv  j  dz 

.  — oo  -»no  0 


x  er  (z)  exp  ( Iku )  exp  ( kz ) 


P(t/  -  /(z-zp)}  -  v(ft  -  /z) 
{«*+v2+(z-z0)^ 


Integration  over  v  then  gives 


?  du  t  dz  ^gfe'igeM  , 

•L  J0  «  +  /<z-z0) 


which  yields. readily  to  contour  integration  over  u  to  give 


(6.38) 


^ri(*o) "  -Mofi’oO  *  X  cosfl  exp(lkuo)  exp  (-te0)  0 

x  J  tr(z)  exp{~2/c(d— z)}  «b  . 

*e 


(6.39) 


For  5^2  we  have 


(6.40) 


x  a(z)  cxpQkH'Tt)  exp(-kz) 


( cos0,sln0,l )  x  (xjr-]f) 
|xo~3f|3 


which  by  a  similar  process  reduces  to 

Bvi’iti'o)  “  k  x  cosS  exp  (/Arw0)  avp  (-Azq)  9 


(6,41) 


x  J  cr(z)  exp{-2/f  (z-zo))  </z  . 

mlnlJoJ) 

For  the  special  case  in  which  the  electrical  conductivity  of  the  seawater  has  the  con¬ 
stant  value  o-o  throughout  the  ocean,  over  the  depth  range  0<zq<cI  the  expressions  (6.39)  and 
(6.41)  for  Byp\  and  fl^>2  respectively  become 


Bmfeo)  “  X  cosfl  exp(-2kd)  exp (iku0)  sink  (kz0)  v 

£w»i(*o)  "  P-& qEou  x  cos0  exp  (- kd )  exp  ( ikua )  sinh {k  (c /— z0))  D  , 


(6.43) 


and  over  the  depth  range  z0>d  they  become 

Bypiix'o)  —  -xh  fjLtfroEoa  cosO  exp  ( iku0 )  exp  ( ~kz0 )  v  (6.44) 

Byp j(*o)  *“  0  •  (6.45) 

For  the  total  field  minus  the  uniform  component  Byy ,  derived  by  adding  together  (6.32) 
and  (6.44),  is  given  for  a  single  term  by 

Sj- Or<j)  "  ^  n&roEoa  exp (iku0)  sxpi-la^ )  sin#  (cos9jln6,i)  . 


(6.46) 
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6.2  Application  to  two-dimensional  bottom  features 

If  we  represent  a  two-dimensional  bottom  feature  in  the  form 

w 

C (a/o)  —  f  a(k)  exp(lku0)  dk  ,  (6.47) 

where  the  a(k)  for  positive  values  of  k  a;e  the  coefficients  for  terms  for  which  — y  <  0  < 

and  the  a(k)  for  negative  values  of  k  are  the  coefficients  for  terms  for  which  0  is  in  the  oppo¬ 
site  direction,  we  then  havs  for  £&  for  zo>0 

Bs  (*o)  “  Vi  Mo^o^o  I  P  j  \fl  ( k)exp  (ikua)+a  (-k)exp  (-  /Aruo)i  exp  (-la o)  dk  (6.48)  . 

o  '  ' 

OP 

+  /  sin0  2  f  la  (k)exp  (ikuQ)-a  (-k)exp  (-ikuQ)]exp  (-1 cz0 )  }  • 

o 

For  0<</<2o  we  have  for  By? \  and  Bypi 

M 

£iwi(*o)  “  ~H(PoEq  COS0  (-sin0,cas0,O)  f  {a  Oc)exp  (/to0)+a  (~k)exp (-iku0))  (6.49) 

of 

exp  (-2kd)sinh  (kz0) 

.  X  1—  exp  (-2kd) 

— 

D^OTo)  "  Mo*? o^o  COS0  (~sin0,coj0,O)  J  {a  (fc)exp  (lku0)+a  (-k)exp  (-Mrw0))  (6.50) 

o 

^  exp  (-fef  )sinh{  A:  (d-z0)}  „ 

1  -exp(-2kd)  ak  ’ 

and  for  z0>d  we  have  for 


£w>(*o)  -  ^triOfo)  -  -'A mo^o^o  cos0  (- si n0, cos 0,0)  J  (6.51) 

*o 

x  (o  (k)exp  (/fcu0)+a  (— A:)exp  (-ikuo)}  exp(—kz§)  dk  . 


The  bottom  features  to  which  this  method  is  most  easily  applied  ate  low  ridges  and 
shallow  trenches,  which  may  be  roughly  approximated  by  a  Gaussian  profile  of  th*;  form 


t  («o)  -  A  exp 


«(? 

2qj 


(6.52) 


ftA  L1*  * 
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where  A  is  positive  for  a  ridge  and  negative  for  a  trough,  and  a  is  the  effective  half-width  of 
the  feature  (see  figure  2.2);  to  properly  satisfy  the  condition  (6,18),  a  should  be  somewhat 
larger  than  A.  If  (6.52)  is  expressed  in  the  form  (6.47),  then  we  have  for  afk) 


a  (k)  -  -4£L.  ex 

V2rr 


f-1 


(6.53) 


When  expression  (6.46)  is  adapted  to  bottom  features  expressed  in  the  form  (6.47),  it  becomes 


•*  {h  fiip-QEaA  {sinfl  (cos0,s/«0,O)  /t(x o/x)  —  sin9  2  /j(xo,a)} 

M 

/iCx^voe)  ~  dT1  f  [a  (k)exp  (iku^)  +  a  (-k)exp  (~/ku0))  exp  ( -kzo )  <Jk 


(6.54) 


(6.55) 


/aOcov «)  -  i  J  {o  (k)exp  ( iku0 )  -  a  (-k)exp  (r-tku0)}  exp(-kzQ)  dk  .  (6.55) 


Substitution  of  (6.53)  gives  for  I\  and  /2 


I\(x9#)  -  ^  expj  -  -S—  -  J  cos(*«y)  dfc 

/2<ac3^K>  «  /  exp  -  -Sy-  -  z<>*  (too)  d*  • 


(6.57) 


(6.58) 


These  integrals  are  dealt  with  as  Fourier  cosine  and  sine  transforms  in  ErdeTyi  et  at.,  1954,  v.  1 
(sec.  1.4,  eq.  16,  p.  15  and  sec.  2.4,  eq.  27,  p.  74),  and  have  the  solutions 


I\(XQfit)  -  V4  !/+Cx0/t)  +  f-  (xq,ck)} 

/20co^a)  —  f/2  (^.Cv^^a)  —  , 

where 

/+  -  exp  (s. I )  erfc  (s*.) 

.  Zq+  too  .  v  , 

5+  *'  ~~yp2a —  ’  '*>0  ’  *0>® 

/._(xo*«)  -  expfj* )  cr/c (s_) 

Z0  “  /tf  o  . 

*-  -  - °  ,  «>0  ,  zo^d  , 

with  the  function  er/cix)  tvi.ng  defined  such  that  on  the  positive  real  axis 


(6.59) 

(6.60) 


(6.61) 


(6.6?.) 


(6.63) 


(6.64) 
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h— r-H 


FIGURE  2.2 


Illustration  of  the  parameters  describing  a  Gaussian  ridge  or  trough  on  the  ocean  floor,  as  used 
in  equation  {6.S2).  A  is  the  ridge  height  (positive  for  a  ridge,  negative  for  a  trough),  a  is  the 
half ‘Width  of  the  rid^e,  (to  is  the  perpendicular  distance  from  the  ridge  axis,  and  {  is  the  vertical 
coordinate  of  the  lidge  surface  a?  a  function  of  tt0. 


(6.65) 


erfcix)  «•  -jmr  J  exp(-t J)  dt  . 

Note  that/+  and/_  are  complex  conjugates  of  each  other,  and  that  therefore  both  (6.59)  and 
(6.60)  will  always  be  real-valued.  The  integrals  /(  and  /2  are  plotted  in  figure  2.3  for  the  special 
cases  of  (a/z0)  equal  to  0.2,  0.4,  and  0.8,  for  values  of  ( u<Ja )  ranging  from  0.0  to  6.0;  I\  is 
symmetric  with  respect  to  uq,  and  /j  is  antisymmetric  with  respect  to  m0.  Expressions  found 
useful  for  the  evaluation  of  (6.61)  and  (6.63)  were  the  formulas  7.1.5  on  page  297  of 
Abramowitz  and  Stegun  and  7.1.23  oh  page  298  of  same,  the  former,  being  appropriate  for  small 
absolute  values  of  the  argument  and  the  latter  for  large  values;  it  was  found  that  if  10  terms 
were  retained  in  the  latter  expansion  then  there  was  a  region  of  overlap  of  satisfactory  size  in 
the  z -plane  where  neither  of  the  corresponding  expressions  for  exp (zJ)er/c(z)  diverged  and 
both  gave  the  same  value  to  several  significant  places. 

As  a  practical  example,  we  consider  the  case  of  a  submarine  canyon  on  the  continental 
shelf.  We  assume  an  electrical  conductivity  for  the  ocean  of  3.3  S/m,  and  as  a  calculational 
convenience  we  consider  the  special  case  in  which  the  feature  is  parallel  to  the  Inducing  EMF 
(such  a  feature  running  perpendicular  to  the  inducing  EMF  yields  a  zero  anomalous  magnetic 
field  above  the  ocean  surface,  although  an  anomalous  field  is  present  below  the  surface).  Also, 
in  order  to  satisfy  the  approximations  made  in  deriving  this  method,  the  feature  must  be  some¬ 
what  broader  than  it  is  deep  and  somewhat  less  in  vertical  dimension  than  the  depth  of  the 
ocean.  We  take  the  value  of  the  Inducing  EMF  to  be  1.0  /iK/m,  a  reasonable  value  for  an 
electromagnetically  quiet  day.  The  submarine  canyon  extends  to  a  depth  20  meters  below  an 
otherwise  level  ocean  floor  at  a  depth  of  100  meters  beneath  the  ocean  surface,  and  has  a  Gaus¬ 
sian  shape  with  a  standard  deviation  of  80  meters.  The  anomalous  magnetic  field  due  to  the 
undersea  canyon  is  plotted  in  figure  2.4,  normalized  against  the  spatially  uniform  magnetic  field 
calculated  from  equation  (2.5)  for  such  an  ocean  and  such  an  inducing  EMF  in  the  absence  of 
any  bottom  feature;  this  field  is  -  2.07x  10“10,p  tesla.  Note  that  while  the  anomalous  field  is,  at 
its  strongest,  no  more  than  about  a  tenth  the  strength  of  the  uniform  field  upon  which  it  is 
superposed,  the  spatial  gradient  properties  of  the  anomalous  field  should  make  it  easy  to  distin¬ 
guish  from  the  uniform  background. 


FIGURE  2.3 


The  integrals  /i  (solid  line)  and  h  (dashed  line)  according  to  equations  (6.57)  and  (6.58),  for 
three  different  ridge  widths.  Notations  are  a  —  ridge  half-width,  d  -  depth  of  ocean,  and>  - 
lateral  distance  from  center  of  ridge. 


The  anomalous  magnetic  held  components  By  and  B,  induced,  at  the  sea  surface  over  a  subma¬ 
rine  canyon  incised  parallel  to  the  x-axis  into  a  level  seafloor,  by  an  electric  held  running  paral¬ 
lel  to  the  axis  of  the  canyon,  as  given  by  equation  (6.53).  The  held  components  are  normalized 
against  the  reference  By  held  which  would  be  induced  by  the  primary  electric  held  in  the 
absence  of  the  canyon,  with  a  completely  level  seafloor.  For  this  example  the  ratio  of  a,  the 
canyon  half-width,  to  d ,  the  depth  on  the  level  part  of  the  ocean  floor,  is  a/(/*-0.8,  and  the 
ratio  of  canyon  depth  to  water  depth  is  0.2.  The  distance  y  is  the  lateral  distance  from  the  axis 
of  the  canyon. 


TyyyryTjrrjrrpT^j*  «>  ,•_* •>  -j.  ■>  r;»  ■>  -.« ~f  -..  r.. 


tiT  v\  V*. v*  ;y*\.  v„  wt.  v\  *r.  w.  *cj  , 
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(.3  Application  to  three-dimensional  bottom  features 

Let  us  now  suppose  that  we  have  a  three-dimensional  bottom  feature  whose  shape  is 
given  by  the  function  {  tx.y),  which  we  write  as  a  Fourier  transform; 


•9  .  m 

{(xj»)  -  J*  dkx  J  dky  aikx,ky)  txviiikxX+kyy)}  . 


(6.66) 


If  the  direction  of  positive  x  is  chosen  to  be  parallel  to  the  uniform  current-inducing  electric 
field  within  the  ocean,  then  we  can  write 


kx-k  cos9  ,  ^  —  Ac  sine  ,  *J-Ar,J.+  fc2  , 


(6.67) 


where  k  is  a  positive  real  scalar,  which  will  transform  the  Cartesian  coordinate  frame  in  which 
(6.66)  is  expressed  into  the  cylindrical  coordinate  system  In  which  (6.46)  is  expressed.  Equa¬ 
tion  (6.46)  is  the  expression  for  the  anomalous  magnetic  field  5J.(*o)  above  the  ocean  surface 
due  to  a  bottom  feature  given,  in  Cartesian  coordinates,  by  a  txp{l  (kxx^-kyyo))\  as  (6.66)  is  a 
superposition  of  such  terms,  we  thus  have  for  ffA  0?)  the  anomalous  field  above  the  ocean  sur¬ 
face  due  to  the  bottom  feature  {(x,y) 


&  0?)  -  J  dkx  J  dky  A  ( kxyky )  txp{t  (kxx+kyy))  exp(-kz > 

•  —99  —99 

it  jz  J ) 
k* 

'  * 

or,  in  cylindrical  coordinates, 

fa  (?)  *  J  k  dk  J*  de  a  (A :  zo$9,k  sinP) 


(6.68) 


(6,69) 


x  exp (ikix  cos®  +  y  sir#)}  exp(-kz)  $in0(co,s0,sin0,/)  . 


In  cases  where  {tx*y)  is  cylindrically  symmetric,  a(k  cos9,k  sinfl)  is  a  function  of  k 
only,  and  judicious  use  of  the  identity 


2 ft  J0(kp)  -  J f  d9  exp(M  (x  cosfl  -(-  y  sin0)j  ,  p 1 «  .v2  +  y1 
o 


(6.70) 


allows  (6.69)  to  be  written  in  the  form 
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/PjOf)  —  -  irju 


{ka(Hexp(  kz)  ^  ^  ^  ^  ®, 


(6.71) 


x  Mkp)  dk  .; 


from  here,  use  of  the  standard  Bessel  function  recurrence  relations  gives 


%  Of)  -  -  ir/i#ro£o  |  k  a  (k)  exp  (■ -te )  I  y2(Ap) ,  £  J2(kp)  -  y,0tpX(v72) 


%Mkp)^dk  . 


If  the  bottom  feature  has  the  Gaussian  form 


£(*,?) -4  eX| 


then  we  have  for  its  Fourier  transform 


f-^l' 


dxt|  — -  -  ~  expj  ~ 


-a(k)  , 


and  (6.72)  becomes 


««  -  -  |  * 

*  I  ^  J,(kp) ,  Mkp)  -  ±  J,«p) ,  j  J,(kp)  ]  dk  , 


which  with  the  substitution 


v~k  p 


turns  into 


where 


I\tX#)  -  f  exp  -  4“  exp  -  —  v  J\(v)dv  , 

T>  P  2  P 


(6.73) 


(6.74) 


(6.75) 


(6,76) 


W__  [  «  /,«„) ,  4  /.of*)  —  /|0?,a)  ,  ^  /2C?,a)  1  ,  (6.77) 

7a*  /\  *  n*  1  n 


(6.78) 
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IiOtfit)  -  ^ 


v  exp - 5- 


a2  v2 


and 


h&fit) "  ^  v 


exp  - 


2  2 


a1  v1 
o1  2 


exp  - 


J i(v)  dv 


(6.79) 


exp  — ■ •* *  v  72(v)  dv 


(6,80) 


These  integrals  may  be  straightforwardly  evaluated  in  the  two  regimes  p  «  a  and 
p  »  a,  In  the  first  case,  the  Gaussian  in  the  integrands  has  gone  to  virtually  zero  while  v>  is 
still  very  small;  accordingly,  if  the  Bessel  functions  in  the  integrands  are  expressed  as  power 
series  expansions,  only  a  small  number  of  terms  from  each  expansion  are  required  to  ade¬ 
quately  approximate  the  values  of  the  integrals.  The  basic  integral  for  these  calculations  is 


S„  Of, a)  «  f  v"  cxpl 
*0 


■»{m!)  &- 
a 


e2  2 


exp 


-  —  v  \  dv 


(6.81) 


hl  ttXO  * 


2 

^2  ^-(ii+i)(t/a) 


(Erde'lyl  et  al„  1954,  v,  l,  sec.  4.5,  eq.  24,  p.  146),  where  n  is  a  nonnegative  Integer  and  the 
D„  Cx)  are  parabolic  cylinder  functions  (see  Erde'lyi  et  al.,  1955,  v.  2,  sec.  8,  p.  115),  and  the 
appropriate  power  series  expansion  for  Bessei  functions  of  integer  order  is 


y',Cx)“*?o  M(fo+*)!|' 


A  »+Jk 
2 


(6.82) 


Accordingly,  from  the  integral  representations  /1,  /2,  and  /j  we  have 


/lOTyt*)  ■“  £  -2*1. 1 


and 


IjOtpt)  -  3C  ^“1 


tiOtv)  -  £ 


*T0  22*4-1  (A !)  { (A-l-l)t) 


■(-!)* 


(A !)  { (A  4- 1)!) 


(-D* 


Sik+\O?0) 


S  2*426? ,«) 


(6.83) 

(6.84) 


Sjk+tOfflt) 


*-0  22*+J  (A !)  { (A4*2) !) 

To  generate  the  necessary  parubolic  cylinder  functions,  we  have  the  recurrence  relation 
/>»+ iCx)  -  x  Dn  0c)  +  n  DH.. |0e)  -  0 


(6.85) 


(6.86) 
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(ibid.,  sec.  8.3,  eq.  14,  p.  119),  and  the  functions 

f)-i(v)  -  -\/^r  exp|~  erfcl-y^)  ,  (6.87) 

and 

0-a<* )  -  axp[-^j  |  1  -  x  exp  ^  trfc  (-fy  (6.88) 

-  exp|-~J  -  x/LjOc)  , 

where  trfc  be),  the  complementary  error  function,  Is  as  given  in  equation  (6.65)  of  the  last  sec¬ 
tion.  Expressions  (6.87)  and  (6.88)  were  found  by  comparison  of  (6.81)  for  the  cases  n» 0  and 
n-1  with  the  equivalent  expressions  in  Grobner  (sec.  312,  eqs.  10(a)  and  10(b),  p.  57). 

From  (6.81)  and  (6.86)  we  have  the  recurrence  relation 
Sn(pri4*)~  (#s—l) ,  (6.89) 

a  tx 

and  from  (6.81)  and  (6.87)  and  (6.88)  respectively  we  have 

-  |^J  Vf  exp(i?)  <6,90) 

•  •  i 

and 

S|<p’,'“)”  t)  1 "  V? (f )  exp(~£t ) erfc  TS-  (6'91) 

•  £  ~  ■Cf  Soipjfit)  . 

a  a* 

In  evaluating  the  complementary  error  function,  the  approximate  expression  (7.1.26)  on  page 
299  of  Abramowltz  and  Stegun,  good  to  high  precision  anywhere  along  the  positive  real  axis, 
was  found  particularly  useful. 

In  the  limit  of  p—0  it  is  seen  that  /j,  /2,  and  /j  all  vanish;  however,  in  the  magnetic 
field  expression  (6.77)  which  uses  these  values  p  appears  in  denominators  taken  to  various 
powers,  and  it  turns  out  that  this  expression  does  not  vanish  in  this  limit.  Referring  to  the 


.vnvvv'  > -:.nv  ,v-  > ■; 
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expansions  (6,83),  (6,84),  and  (6,85),  we  insert  them  into  (6,77)  and  discard  any  resulting 
terms  in  positive  powers  of  p,  to  get 


nacoa— isSaogt 


p*~0 


Slipup) ,  S yip ^2  fit)  -  t  (6,92) 

4p*  4p*  pi 


-4s2(p^^) 

P 


where  the  indicated  powers  of  p  in  the  denominators  cancel  the  p-dependences  of  the  accom¬ 
panying  S„  (sec  (6.81));  since  x  andy  will  vanish  withp,  this  reduces  further  to 


-Vfli 


exp  -p 


er/c\ 


7£l 


(6.93) 


At  the  other  extreme,  where  p  »  a,  the  value  of  the  Gaussian  is  still  close  to  unity 
by  the  time  that  v  is  large  enough  to  make  the  Bessel  functions  for  practical  purposes  vanish;  in 
this  case  we  expand  the  Gaussian  in  a  power  series  about  zero,  and  integrate  termwise.  The 
basic  integral  is 


4(v)  cfv 


I 


which  for  n  and  m  both  nonnegative  integers  and  n  >  m  is 


(6.94) 


(-t)m 


(«  -  w)!  pm  1± 
/3"+l  "  hs  p 


(6.95) 


/3 


1+' 


V2 


(6.96) 


(ErdeTyi  et  al.,  1954,  v,  2,  see.  8.6,  eq.  6,  p.  29  and  Erde'lyi  et  at.,  1955,  v.  1,  sec.  3,4,  cq,  17, 
p.144),  where  the  Pf(x)  are  Legendre  functions  of  integer  degree  n  and  integer  order  m  (ibid., 
sec.  3.4,  p.  143).  The  Gaussian  series  expansion  is 


y  (-1)" 

„To  T  n ! 


a_ 

P 


2h  Yin 
•*  » 


(6.97) 


giving  for  /|,  /],  and  /j 


lllfJyZflt)  - 


f  (-!)« 

*-0  2*  (,1c  I) 


\2k 

~  L} A.lW)  * 

P 


(6.98) 
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Ilip/P.)  -  52  ^  £»+u(p«*)  »  <6.99) 

and 

v2* 

lifpjfit)  “  £  2*(^!)'  ^2*+t.2<P^)  •  (6.100) 

Expressions  for  lo.i  and  1 1,2  are  required  by  (6.98)  and  (6.100)  respectively,  but  are  not  given 
by  (6.95);  from  integral  tables,  these  are 

I0,i(p^)-  1-4-  (6.101) 

PP 

(Erde'lyi  et  al'.,  1954,  v.  2,  sec.  8.4,  eq.  6,  p.  19)  and 
I '  ,  W  ,  \i 

2+^  1-^1  (6.102) 

(Erde'lyi  et  al.,  1954,  v.  1,  sec.  4.14,  eq,  2,  p,  182),  where /8  is  as  given  by  (6.96), 

Useful  recurrence  relations  and  bootstrap  functions  for  computing  the  Legendre  func-  •  * 
tions  required  by  (6.95)  are 

Oi-m+l)fft,  0c)  -  (2«+l)x /»,?(*)  -  («+m)P*.,  Cx)  (6.103) 

and 


(l-“Xa)vaJ^7,+1(x)  —  (n‘-m)xPH> Cx)-  (n+m)  PZL\(x)  (6.104) 


(Erde'lyi  et  al.,  1955,  v,  1,  sec,  3,8,  cqs.  12  and  17,  p.  161),  and 

PS  (x)-l  (6.105) 
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Pi  Oc)  -  3  (l-jc2)  .  (6.110) 

Reference  to  (6.95),  (6.103),  and  (6.104)  gives  the  recurrence  relations  for  L„>m  (p,z) 

(2/if-l)  b/p)L„-U\(pj)  -  n(n-2) 4_2ti  (6.111) 


s 


m 


(n+ 1) \<p ^)  -  (z/p)IHli(p^)  , 


(6.112) 


and  some  useful  bootstrap  functions,  from  (6.95)  and  the  previously  listed  Legendre  functions, 


^1,1  (pi*)  “  pj 


(6.113) 


(6.114) 


note  that  the  recurrcnca  relations  (6.111)  and  (6.112)  are  invalid  with  the  use  of  bootstrap 
functions  for  which  m>tt,  making  it  Improper  to  use  either  (6.101)  or  (6.102)  to  start  them 
off. 


In  the  limit  of  the  integrals  (6.78),  (6.79),  and  (6.80)  reduce  to 

m 

/|“  J  J\(v)dv  -  1  , 


(6.115) 


/j  -  J  vjx(v)dv  «  1  , 


(6.116) 


•o 

/j-  J  vJ2(v)dv  -  2  ; 


(6.117) 


these  integrals  are  special  cases  of  /.o,i,  Li,i,  and  L  1,2,  and  their  values  were  derived  from  the 
appropriate  expressions  in  the  text  (expressions  (6.101),  (6.113),  and  (6.102)  respectively)  by 
taking  the  appropriate  limits,  Substitution  of  these  values  into  (6.77)  gives 


lim  Ba  (3?) 

p.  -*o 


_  <*  1  7xy  2x1  _  1  2. 

m  «  3  '  3  >  ~  7  ~  > 


P2  ’  P2 


(6.118) 


Computation  of  the  elements  of  the  tensor  gradient  of  the  anomalous  magnetic  field  is 
.  expedited  by  the  fact  that  spatial  derivatives  .of  the  members  of  the  sets  of  the  S„  and  the  Ln  m 
are  also  members  of  the  respective  sets;  from  the  integral  expressions  (6.81)  and  (6.94)  we 
have 

-  -^Sn+\(pjjz)  +  2jS„+2(pjfit)  ,  (6.119) 

Op  p*  p 

j^S„<p,Zfic)  -  -  — Sa+itp^a)  ,  (6.120) 

8  r 

-  —  ^rt+t,*(p/.a)  .  (6.121) 

and 

(p£fx)  ™  (p^pi)  .  (6.122) 

Modification  of  the  routines  used  to  calculate  I\,  /2,  and  /3  for  the  calculation  of  their  spatial 
derivatives  should  be  fairly  straightforward. 

A  prabtical  problem  in  the  use  of  the  wcutr-'t'cc  relations  (6,89)  and  (6.111)  in  the 
calculation  of  I\,  I2,  and  13  is  that  for  many  choices  of  p,  r,  and  a  the  series  diverge, 
apparently  ta  a  result  of  numerical  instability;  where  this  occurred,  investigation  of  thr  problem 
showed  that  successive  terms  in  the  expansions  for  7i,  /j,  and  I3  first  decreased  in  magnitude 
to  values  very  small  in  comparison  to  the  magnitude  of  the  first  term  in  the  series,  and  then 
increased  without  bound.  It  was  found  that  if  the  expansions  were  truncated  just  before  the 
point  of  divergence  was  reached,  the  agreement  between  the  equivalent  series  in  the  S„  and  the 
Lnjn  in  their  areas  of  mutual  validity  was  good  to  several  significant  places,  suggesting  that  this 
procedure  is  a  workable  solution  to  the  divergence  problem.  However,  in  many  cases  it  will  be 
desirable  to  get  at  least  a  rough  estimate  of  the  number  of  terms  actually  required  for  a  result 
of  a  given  level  of  accuracy,  either  to  avoid  computing  unnecessary  terms  or  to  determine 
whether  or  not  a  series  expansion  is  invalid  due  to  divergence  before  the  requisite  number  is 
reached.  Both  methods  are  based  on  the  series  expansion  in  ascending  powers  of  its  argument 
of  a  factor  in  the  integrand  of  an  integral,  the  justification  being  that  another  factor  in  the 
integrand  will  have  effectively  vanished  before  the  argument  of  the  expanded  function  becomes 
large  enough  to  invalidate  its  series  expansion.  To  produce  a  series  truncation  error  estimate 
for  a  given  choice  of  p,  z,  and  a,  one  first  truncates  the  integration  at  a  selected  cutoff  point, 
and  estimates  the  error  introduced  by  disregarding  the  remainder  of  the  integral;  then,  one  cal¬ 
culates  the  difference  at  this  cutoff  point  between  the  value  of  the  factor  tr>  be  expanded  and 
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the  value  of  its  expansion,  truncated  after  the  specified  number  of  terms.  A  reasonably  reliable 
upper  limit  on  the  error  produced  by  the  series  truncation  may  be  arrived  at  by  making  the  con¬ 
servative  assumption  that  the  percentage  discrepancy  between  the  factor  value  and  the  value  of 
its  expansion  is  the  same  over  the  entire  range  of  integration  from  zero  up  to  the  cutoff  point, 
and  the  total  error  is  calculated  by  multiplying  this  percentage  error  by  the  value  of  the  trun¬ 
cated  integral  based  on  the  truncated  series  expansion  and  adding  the  result  to  the  estimated 
error  from  the  truncation  of  the  integral,  This  process  involves  some  tradeoff  decisions,  in  that 
selecting  a  greater  cutoff  point  will  decrease  the  integration  truncation  error  but  increase  the 
scries  truncation  error  unless  a  larger  number  of  terms  is  used  in  the  expansion. 

As  an  example,  we  consider  the  case  of  a  Gaussian  hill  20  meters  high  and  of  width 
such  that  a  ■*  80  meters,  rising  from  an  otherwise  level  ocean  floor  100  meters  below  the  sur¬ 
face.  The  magnetic  field  is  computed  at  the  points  of  a  5x  5  horizontal  grid  at  the  ocean  surface, 
centered  directly  over  the  peak  of  the  hill  and  with  a  grid  spacing  of  160  meters,  with  faces 
parallel  to  the  direction  of  the  inducing  EMF.  The  water  depth  is  not  untypical  of  a  continental 
shelf,  and  underwater  features  such  as  a  reef,  a  dune,  or  a  highly  eroded  suomerged  island 
might  have  dimensions  on  the  order  of  this  Gaussian  hill.  A  uniform  inducing  electric  field  of 
1.0  n  V/m  in  the  £  direction  is  assumed,  and  the  ocean's  electrical  conductivity  is  taken  to  be  a 
uniform  3.3  S/m.  Note  that  most  of  the  points  on  the  grid  are  in  the  intermediate  region 
between  p  «  a  and  p  »  a\  in  this  region  both  of  the  series  expansion  schemes  given  above 
should  be  usable,  and  in  fact  for  this  problem  there  was  a  narrow  range  of  values  of  p  where 
they  gave  the  same  values  for  /j,  /21  and  1 3  to  several  significant  places.  However,  difficulties 
encountered  suggest  that,  as  a  safety  precaution,  for  computations  with  any  given  combination 
of  «  and  z  one  should  compare  the  outputs  of  tl  ‘>  two  algorithms  over  the  range  of  p  with 
which  one  is  concerned,  in  order  to  establish  whether  there  is  in  fact  a  range  of  equivalent  out¬ 
put,  and,  if  so,  exactly  where  it  is.  The  results  of  the  calculations  are  presented  in  table  2.3. 

The  scaling  rules  discussed  toward  the  end  of  section  3  may  be  used  with  these  results 
to  derive  the  magnetic  values  at  the  ocean  surface  for  any  situation  where  the  relative  dimen¬ 
sions  of  ocean  depth,  hill  height,  hill  width,  and  grid  spacing  are  the  same  as  in  this  example: 
Also,  the  field  values  are  directly  proportional  to  both  the  strength  of  the  inducing  EMF  and 
the  electrical  conductivity  of  the  seawater.  According  to  equation  (6.77),  the  field  values  are 
directly  proportional  to  the  height  of  the  hill,  but  it  should  be  recalled  that  this  equation  loses 
its  validity  if  the  height  of  the  hill  approaches  either  half  of  the  ocean  depth  or  the  width  of  the 
hill. 
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X 

-320.0  m 

.  -160.0  m 

0.0  m 

160.0  m 

320,0  m 

/- 320.0m 

Bx 

0.8S0 

0.926 

0.000 

‘  -0.926 

-0.850 

B, 

0.153 

-0.377 

-0.851 

-0.377 

0.153 

2 

-0.884 

-1.732 

-2.357 

-1.732 

*0;884 

/-  160.0  m 

Bx 

0.926 

1.543 

0.000 

-1.548 

-0.926 

By 

1.011 

■  1.205 

1.044 

1.205 

1.011 

B 

-0.866 

-2.666 

-4.720 

-2,666 

-0.866 

/-0.0m 

Bx 

.0.000 

0.000 

0.000 

0.000 

0.000 

By  . 

1.748 

3.788 

9.569 

3.788 

1.748 

B, 

0.000 

0.000 

0.000 

0.000 

O.COO 

/--160.0m 

Bx 

-0.926 

-1.548 

0.000, 

1.548 

0.926 

By 

1.011 

1.205 

1,044 

1.205 

1.011 

B, 

0.866 

2.666 

4.720 

2.666 

0.866 

/——320.0  m 

Bx 

•0,850 

-0.926 

0.000 

0.926 

0.850 

By 

0.153 

-0.377 

-0.851 

-0.377 

0.153 

Bx 

0.884 

1.732 

2.357 

1.732 

0.884 

TABLE  2.3 

Held  values  at  the  ocean  surface  over  a  submerged  Gaussian  hill  as  described  in  the  text  on 
page  67.  The  field  values  are  in  units  of  \(Tn tesla,  and  the  origin  of  the  coordinate  system  for 
position  on  the  surface  is  on  the  vertical  axis  of  symmetry  of  the  hill.  In  the  absence  of  the  hill 
the  uniform  field  at  and  above  the  surface  would  be  B — 2.07x  1(T10^  tesla . 


7.  Interucttaa  with  the  Earth’s  magnetic  held  of  seawater  current  flow 

An  electric  current  will  be  generated  by  the  flow  of  seawater  through  the  Earth’s  mag¬ 
netic  field,  and  this  electrical  current  will  in  turn  generate  a  magnetic  field  of  its  own.  As  well 
as  the  magnetic  anomalies  associated  with  the  flow  of  tides  and  steady  ocean  currents  and  with 
turbulent  water  motion,  Which  are  outside  of  the  scope  of  this  study,  there  will  also  be 
anomalies  Mspci^tud  with  the  diversion  of  water  flow  over  the  ocean  bottom  by  topographical 
features  of  the  ocean  floor.  Unfortunately,  the  flow  patterns  of  real,  viscous  water  about  even 
very  simple  shapes  are  difficult  to  model,  and  the  subject  has  in  general  not  been  well  studied; 
accoidlngly,  even  though  this  effect  is  probably  significant  in  comparison  with  the  others  dis¬ 
cussed  in  this  report,  it  cannot  be  covered  here  because  the  available  data  la  insufficient  even  to 
estimate  orders  of  magnitude.  For  some  examples  see  Cox  and  Sandstrom,  1962,  and  Osborne, 
1980. 
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Chapter  III 


Magnetic  Anomalies  Induced  by  Fronts  within  the  Ocean 

1.  Introduction 

Oceanic  fronts  are  regions  where  dissimilar  water  types  abut  one  another,  and  they  can 
exist  as. stable  conditions  because  the  electrical  conductivity  of  seawater  increases  both  with 
increasing  temperature  and  increasing  salinity,  whereas  the  density  of  seawater  increases  with 
increasing  salinity  but  decreasing  temperature.  As  a  result  of  this,  two  bodies  of  water  can  lie 
adjoining  each  other  without  a  density  contrast  but  with  appreciable  contrast  in  electrical  con* 
ductivity,  and  fronts  with  their  main  contrast  in  the  vertical  niay  have  density  as  well  as  con¬ 
ductivity  jumps.  The  causes  of  temperature  and  salinity  variations  lie  ultimately  in  air-sea 
processes  such  as  precipitation,  evaporation,  insolation,  and  long  wave  back  radiation,  but  the 
actual  distribution  of  water  types  within  the  ocean  results  from  dynamic  processes  which  tran¬ 
sport  regions  of  water  that  have  been  subjected  to  these  influences  in  different  proportions. 

Some  of  the  more  common  of  these  dynamical  processes  derive  from*  gravity  waves.  In 
this  instance,  we  regard  the  redistribution  of  water  types  in  these  so-called  internal  waves  as  a 
form  of  "frcntaT  creation*,  these  fronts  have  time  scales  in  the  range  from  the  inertial  period  to 
the  Vaisala  period  (Eckart,  1960).  Other  processes  -  currents  and  turbulence  -  can  lead  to 
longer  time  scales.  Internal  waves  are  ubiquitous  in  the  ocean  but  have  local  intensification  in 
shallow  water  (Oordon,  1977).  Since  the  vertical  structure  of  the  ocean  is  always  one  of 
increasing  density  with  depth  (usually  caused  by  decreasing  temperature),  internal  wave  fronts 
are  always  present.  The  more  persistent  fronts  where  density  remains  almost  constant  demand 
the  close  proximity  of  contrasting  salinity/temperature  water  types. 

This  chapter  deals  with  the  calculation  of  the  magnetic  anomalies  associated  with  the 
electrical  conductivity  properties  of  oceanic  fronts.  Section  2  uses  a  perturbation  approach  to 
derive  equations  describing  the  alteration  of  the  electric  field  within  the  ocean  by  the  presence 
of  a  front,  and  discusses  the  calculation  of  the  magnetic  effects  of  the  front  by  use  of  the  solu¬ 
tions  to  these  equations.  Section  3  deals  with  solutions  to  the  electric  field  equations  for  a  spe¬ 
cial  case  of  particular  interest,  section  4  illustrates  the  practical  use  of  the  work  of  section  3, 
and  section  S  gives  two  different  methods  for  a  general  solution  for  the  electric  fields  in  the 
perturbation  limit,  one  being  the  more  convenient  for  analysis  of  Fourier  transformed  data  and 
the  other  being  the  more  suitable  for  the  study  of  localized  fronts  in  real  space. 
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2.  Basic  equations 

2.1  The  oceanic  electric  field 

We  start  by  supposing  that  the  front  exists  as  a  steady  state  condition,  making  the 
magnetic  field  associated  with  it  constant  in  time.  From  Maxwell's  equations  we  then  have  for 


£  the  electric  field  within  the  ocean 

Vx£*-0  ,  (2.1) 

which  implies  that  £  may  be  expressed  as  the  gradient  of  some  scalar  function 

£-V*  .  (2.2)  • 

For  7  the  electrical  current  density  distribution  within  the  ocean  we  have 

7  -«r£  ,  (2.3) 

where  <r  is  the  electrical  conductivity  distribution  of  the  seawater,  and 

V-7-0  ,  (2.4) 

another  consequence  of  the  steady  state  assumption.  Combining  (2.2),  (2.3),  and  (2.4)  gives 
V-<*rV*)-0  .  ,  (2.5) 

or 

or  Vty  +  (V<r)(V*)  -  0  .  (2.6) 


For  fronts  in  the  ocean,  the  variations  from  the  mean  electrical  conductivity  of  the  frontal 
region  are  seldom  more  than  a  few  percent  of  the  mean  conductivity,  suggesting  the  use  of  per- 
.  turbation  methods  in  solving  (2.6).  Accordingly,  we  write  or,  which  we  may  take  as  a  given,  in 
the  form 

trCf)  **  <To(r)  4-  *  crjC?)  (2.7) 

and  <j>  in  the  form 

d»Of)~^oC5?)  +  «^iO?)  +  «J^Of)4- •••  ,  (2.8) 

where  e  is  a  perturbation  parameter  that  will  eventually  be  set  to  unity,  and  crt  is  a  small  per¬ 
turbation  on  the  mean  conductivity  cr0,  which  is  taken  for  mathematical  convenience  to  be  hor¬ 
izontally  stratified.  Substitution  of  (2.7)  and  (2.8)  into  (2.6)  gives 
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<roVtyo+  (VcroVCV^o) 

+  «  (o"o  +•  <^1  VJ^o  +  (V«ro)-(Vd>t)  +  (Vo-i)*(V^o)) 

+  «J  {^o  Vtyj  +  «rj  V2^i  +  (VtroJ'CV^j)  +  (Vo,i)*(V^i)J 
.+  •••  -0 .  .  * 

Requiring  that  (2.9)  hold  for  till  choices  of  the  value  of «  yields 
<r0V^ o+ (Vo-o)  C7^0)-0  , 
cro  VVt  +  (VoroJ’C^i)  +  <r\  V2^o+  (V«ri)‘(V^o)  “  0  ■, 
<r0V^2+  (Vo'o)‘(V^j)  +  ar\  Vfyj  -f  (VcrjVCV^i)  ■»  0  , 
and,  for  higher  orders, 

0-oV^A  +  (Vcro)-(V0#)  +  (r,  (Wi)'CV<*n_t)-0  . 

Making  the  definition 


/(*)- 


1 

<ro(iiT 


dao(z) 


gives  from  this 


VVo  +  /CzH-V*o-0  , 

V^ki  +  f(z)  2-V0J  -  pjCSf) 

where 

.  ^  {or,0?)/(z)2-V<r1C?))-V*o 

p,a) - - - 

V^2  +  /(z)  2-V^j-^ Of) 


where 


{<riC?)/(Z)2- V<t,0?)}  •  V0,  crjO?)  ^ 

',m - ^U) - 


and 

V^+ZW-V^-^C?) 


(2.9) 


(2.10) 

(2.11) 

(2.12) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 

(2.19) , 

(2.20) 
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*w - -^U) - ^  ',-,c?)  ;•  <ul' 

Setting  c  to  1,  we  have 

<rOt)  —  ctq(z)  •+  cjC?)  (2.22) 

and 

d(3f)-^oO?)  +  ^iO?)  +  ^2(S')+  *  .  •  (2*23) 

As  ld]l  will  be  smaller  than  |0tl  by  something  bn  the  order  of  magnitude  of  the  ratio  of 
(max  |cr  j  |}  to  |tr0|,  it  will  generally  be  practical  to  ignore  all  terms  beyond  <f>\  in  (2.23). 

As  an  important  special  case,  we  suppose  that 

0oC?)-£ox  ,  (2.24) 

which  is  seen  to  be  a  valid  solution  to  (2.1$);  then,  since  0—0O  if  is  uniformly  equal  to 
zero,  in  the  absence  of  a  front  the  electric  field  2?  throughout  the  ocean  is 

'E-V<t>0"E0&  ,  (2.25) 

which  we  may  take  to  be  induced  by  some  outside  source  of  EMF.  If  trt  is  nonzero,  then  we 


have  for  from  (2.16)  and  (2.17) 

. 

V^iCf)  +  /<*>  -  Pifr) 

(2.26) 

_  £0  W 5?) 

»x  • 

(2.27) 

The  choice  (2.24)  for  0C  is  appropriate  for  virtually  all  cases  of  practical  interest  dealing  with 
fronts  near  the  surface  of  an  ocean. 
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2.2  The  Fourier  transformed  magnetic  field 

The  Biot-Savart  expression  for  the  magnetic  field  f?  due  to  a  steady  state  electrical 
current  density  field  7  is 


if  we  define 


PC?)- 


Me  3? 


4»r  js?  I3  ’  ' 

we  note  that  (2.28)  becomes  a  convolution  integral  of  the  form 

B(^)-  f  7<SC)xP<%r*)<*  • 

Y 

Defining  the  Fourier  transforms  S(Jc ),£(£),  and  3  (£), 

7c?>  “/*<£)«**<£  , 

and 

P<X)-f  &<&'**&  , 

we  find  that  (2.30)  is  equivalent  to 

?(D  -  (2tr)3  RQc)x  3(£)  . 

The  current  density  field  7  has  the  form 
7C?)-<rO?)  V*0?)  ; 
if  we  make  the  definitions 

<rC ?)- J*X(E)  e®*  die 
and 

<t>(X)  «•  J  4>(£)  e1*'*  &  , 

we  find  that  (2.35)  is  equivalent  to 

J?(F)-J  iP2$~i?)'b(k').dP  . 


(2.28) 


(2.29) 


(2.30) 


(2.30 

(2.32) 

(2.33) 

(2.34) 

(2.35) 

(2.36) 

(2.37) 

(2.38) 


(2.39) 


i 

i 


i 

i 

i 

i 
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An  equivalent  expression  to  (2.29)  for  P(X)  is 


PC?)~- 


MO 

4jt 


some  work  with  this  expression  shows  that ... 

<2*)-3  R  u$)  ,  (2.40) 

where 

irOr) -/!«*•**<*  , 
giving  for  (2.34) 

where 

T  (It)  ~  l  Tex  ROD  (2.43) 

-  J  (Px  Tc)  ZOc-P)  Q(P)  dP  . 

If  expression  (2.41)  is  written  in  Cartesian  coordinates,  then  use  of  symmetry  about  the  origin 
reduces  it  to  three  sequential  cosine  transforms,  which  are  given  in  Erde'iyi  et  at.,  1954,  v.  1 
(sec.  1.3,  eq.  7,  p.  11,  sec.  1.13,  eq.  43,  p.  56,  and  sec.  1.4,  eq.  1,  p.  14);  the  result  is 

U(j£)-&  .  (2.44) 

k‘k 

The  current  density  distribution  J\  up  to  the  first  order  is  given  by 
j;OP)-ao(z)  Vtfo0?)  +  4  ,  (2.45) 

where 

/JOrt-tnO?)  V*00?)  +  <70(z)  V«^iO?)  ;  (2.46) 

the  first  term  of  (2.45)  is  unaffected  if  <r  j  is  set  uniformly  equal  to  zero,  whereas  vanishes, 
thereby  identifying  itself  as  the  anomalous  part  of  the  distribution,  The  Fourier  transform 
£(£)  of  the  magnetic  field  resulting  from  /<  C?)  is 


(2.41) 


(2.42) 


?(*)-—-  (/(£)  c  (£) 


CA  Qt)  -  J*  (i?x  ic)  {Et(£-F') $<)(/?)  +  io<S-Do,(P))  dP  , 


(2.48) 


where 


<r0(r)  -  J*  Eo(^)  dk  , 

o*i(?)  -  J  I|(E)  , 

^oO?)  “  J*  4>o(£)  $  » 

/  4»i(E)  JP  . 


(2.49) 


(2.50) 


(2.51) 


(2.52) 


7? 


3.  Exponential  decline  of  the  lateral  mean  electrical  conductivity  with  depth 

Let  us  suppose  that  tr0 (z)  the  unperturbed  electrical  conductivity  distribution  of  the 


ocean  goes  as 


1  rf<roU) 
*/oGO  dz 


-  X  ,  X  >  0  , 


where  X  is  a  constant,  giving  an  exponential  decline  in  <r0  with  depth.  Then  (2. IS)  becomes 


VVo  +  X^-0  , 

for  which  (2.24)  is  still  a  solution,  and  (2.16)  and  (2,17)  become 
Vtyt  +  X  -  P\  i- 


where 


Pi  Of)  ■ 


{<r|.X  2  —  V<ri)  •  V^o 


within  the  ocean  and  Is  undefined  elsewhere.  In  this  special  case,  if<60  and  are  provided, 
then  an  analytical  solution  to  (3.3)  Is  possible. 

Going  to  the  Fourier  transform  domain,  we  define 

-  J  Q\0c)  t®*  df<  (3.5) 


PiC?)  -  JT  P\(ft)  dft\  (3.6) 

to  avoid  ambiguity  in  the  Fourier  transform  of  piOf),  we  sot  it  to  zero  outside  of  the  ocean. 
Making  these  substitutions  into  (3  3)  gives 


0(£)<ME)-  PiOt)  , 


where 

Q$)-~V-V+  l\Z-5  . 

From  (3.7)  we  have 


Wv\' 
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whose  inverse  Fourier  transform  4m/>(7)  Is  a  particular  solution  to  equation  (3.3).  Equation 
(3.3)  with  pi  set  uniformly  equal  to  zero  is  solved  by  any  term  of  the  form  «?*'*  for  which  k  is 


such  that  (?(£)- 0.  If  we  write  k  in  the  form 

£-(**.*„/*>  ,  (3.10) 

then  (3.8)  becomes 

Q(Jc)-  k2-kk-k$  ,  (3.11) 

where 

+*,:  1  (3.12) 

setting  Q(Jc)  uniformly  equal  to  zero  gives  for  k 

£  +  k$  .  (3.13) 

Accordingly,  we  have  for  the  general  solution  to  equation  (3.3) 

*,0rWi,<*)  (3.14) 


+  J*  /  A(kx\ky)  exdt(kxx+kyy)-kAz]  dkx  dky 

— m 

to  *o 

+  jf  J  B(kx>ky)  exp {l(kxx+kyy)-'lcBz)  dkx  uky  > 

where  kA  is  the  larger  of  the  two  roots  (3.13)  and  kg  is  the  smaller;  the  integral  term  in  A  is 
downward-increasing,  whereas  the  integral  term  in  B  is  downward-decreasing. 

If  we  suppose  that  A  has  the  same  value  everywhere  from  the  ocean  floor  to  the  sur¬ 
face,  then  the  boundary  conditions  on  any  solution  to  (3.3)  are  that  2'Vt£i-0  both  at  the  ocean 
bottom  and  the  ocean  surface,  which  is  equivalent  to  saying  that  there  is  no  flow  of  electrical 
current  through  either  bottom  or  surface,  This  approximation  is  usually  valid  because  the  con¬ 
ductivity  of  the  seafloor,  while  variable,  is  usually  far  smaller  than  that  of  the  seawater  above. 
From  (3.14)  we  have  that 

to  m  * 

2-V^|CS()  -  J*  f  V(kx,ky>z)  6Kp[i(kxx+kyy)}  dkx  dky  ,  (3.15) 

-to  —to 

VikxJfy*)  “  R  (kx,kyj)  -  kAe~kAtA  (*,,*,,)  -  kBe~ks‘B(k,,ky)  ,  (3.16) 
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1 


4a 


i  ■ 


where 


m 

R(kx>kyj)<-  f  ik,  <t>lP(k)  eV  dk,  . 

—m 

The  boundary  conditions  at  surface  and  bottom  are  seen  to  be  equivalent  to 

,*,*)-  K</Wfriza.)-0  , 


(3.17) 


(3.18) 


where  z5  and  zB  are  the  altitudes  of  the  surface  and  the  bottom  respectively.  These  represent 
two  equations  in  A  and  B,  which  may  be  solved  simultaneously  for  A  and  B  in  terms  of  R  at 
the  two  altitudes, 

In  some  cases  it  may  prove  convenient  to  minimize  the  maximum  value  of  |o-il  by 
constructing  <r0  from  two  different  exponential  curved  joined  together  at  some  depth  within  the 
ocean.  If  this  is  dune,  then  then  the  conditions  (3.18)  still  apply,  and  there  are  additional  condi¬ 
tions  at  the  interfacing  depth,  derived  from  the  requirement  that  the  electrical  current  density 
distribution  be  continuous  across  the  interface;  provided  that  <r o  is  continuous  across  the  inter¬ 
face,  this  is  equivalent  to  saying  that  Vd>i  is  continuous  across  the  interface.  We  already  have 
in  equation  (3.15)  an  expression  for  the  vertical  component  of  the  gradient;  from  equation 
(3.14)  we  have  for  the  horizontal  components 


m  to 

VtfduC?)-  J  J  Kk^kyfi)  H(kit,kyt)  exp(/(^.x+fcy>'))  dkx  dky  , 

H(kx,fcy,z)  “  S(kx,kyj)  +  fkA'A  (kx,ky)  +  e~kaIBQcxtky)  , 


(3.19) 


(3.20) 


where 


m 

SOcxtkyj)-  J  d>i pdt)  elk,‘  dk,  , 


Let  the  lower  region  be  called  region  1,  and  the  upper  regie  n  be  called  region  2,  with  their 
interface  being  at  the  altitude  z/;  bur  object  Is  to  find  in  the  two  regions  by  solving  for  A  i 
and  B\  in  the  lower  region  and  A  j  and  Bi  in  the  upper  region.  For  boundary  conditions  we 
have 


K,(*„*„z*)-0  , 

V\(kx>ky*)-  y2(kx,ky,zj)  , 
// ]  (kx  ^ky  \Z/ )  ■*  kf^(kx  kv  ^Z\ )  , 


(3.22) 


(3.23) 


(3.24) 


S 
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Vi(kxtkyjs)-  0  ,  (3.25) 

giving  one  equation  in  A\  and  B one  equation  in  A 2  and  Bi,  and  two  equations  in  all  four 
parameters.  Although  the  process  is  cumbersome,  theso  conditions  are  sufficient  for  the  unique 
simultaneous  solutions  for  A  j ,  B  \ ,  A  2,  and  B2. 
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4.  Practical  use  of  the  method  of  section  3 

The  basic  approach  used  here  will  be  to  first  compute  the  magnetic  anomaly  field  asso¬ 
ciated  with  a  unit  point  variation  front  the  mean  electrical  conductivity  distribution  a-0,  and  then 
use  this  as  a  Green’s  function  to  get  the  field  for  an  arbitrary  variation  distribution.  Drawing 
upon  section  2.2,  we  have  that  if  By-  is  the  magnetic  field  associated  with  a  variation  distribu¬ 
tion  or  i,  then  its  Fourier  transform  Cy  is  given  by 

5(tl-*S»CIEl  ,  (4.1) 

UG)--$gg  ,  (4.2). 

kk 

EyQi)  •*  i  kx  KyQc)  ,  (4.3) 

where  Ky  is  the  Fourier  transform  of  Jy,  the  anomalous  electrical  current  density  distribution 
associated  with  the  variation  distribution; 


/p0?)-/  /rv(k)  e^dk  . 
From  equation  (2.46).  we  have  that 


Jy(%)  ■■  /kiC?)  +  » 


where 


JviC?)  -  *,(?)  v^oC?) 


Jvi&)  -  o-0(?)  V<£iG?)  ,  (4.7) 

with  the  Oth  and  hr  order  electric  field  potentials  <f>0  and  4>\  as  defined  in  section  2.1.  We 
assume  that  4>0  is  given  by 

*o0?)-£o*  .  (4.8) 

so  as  to  give  a  uniform  electric  field  of  £0  *  within  the  ocean  in  the  absence  of  a  variation  dis¬ 
tribution,  and  that  the  mean  distribution  cr0  is  given  by 

oroOf)  —  W(z)  0-0(24)  exp(\(zs-r)]  ,  (4.9) 


1  ,  zB  <  1  <  zs 
Wv)  -  0  ,  otherwise  , 


(4.10) 
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where  X  is  a  constant  greater  than  or  equal  to  zero,  zs  is  the  altitude  of  the  surface  of  the 
ocean,  and  za  is  the  altitude  of  the  ocean  bottom.  Assuming  for  the  moment  that 


t'lC?)-  8(?-io)  , 

where  8  is  the  Dirac  delta  distribution,  written  as 


m  «•  m 

8(*) -  -srrr'J  dkx  /  dk,  J*  dk,  exp  (-/£:?)  , 


(2ir)J  L 
we  immediately  have  for  Afpi 


JfCyiOc)  «■  X 


i?0 


(2 it)3 

and  hence 

Qcxi) 

For  fTvi  we  have  that 


exp  (-Sc'Xo)  , 


Eo 


(2ir)3 


exp(-iic-)T0)  . 


KrzOc)  “  ~^jr  J  dx  ] '  <fy  J dz  exp(-licX) 


x  W'(*)o-o(zs)exp{X(zs-z)}  V^jO?)  ; 


if  is  the  Fourier  transform  of  0jOf),  then  this  becomes 

<,0fe)«r(-xzs)  J 


Kvidc)  “ 


2ir 

exp {t(k,'-k,-i\)zs)  -  exp [l(k,'-kz-i\)zB} 


kf'-kj-iX 


from  which  we  have 


Cvl(k)  -  , 


.  cr0(2T)exp(-Xzs)  tT 


2n 


(kxi)  J  a,-v)4>i(A*A  ,v) 


e\p{l(k,'-kI-iX)zs)  ~  exp{ /  (Ar  '~ks~l\)zB\ 


kt'-k,-ik 


(4.11) 


(4.12) 


(4.13) 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


To  determine  Oj,  we  refer  to  section  3.  From  equation  (3.14)  we  have  that 
Ql(ic)-*xF(£)  +  A(k?!,ky)6(kI-tkA)  +  B(kXiky).8(kt-ikB)  ,  (4.18) 


where 


kA  -  a  +  0  , 
kB  -  a  -  0  , 


-*• 


(4.19) 

(4.20) 

(4.21) 


and 


P 


-yPx+Oi 


-0->O  , 


where 

ki-k}+ky2  . 

(?)  is  given  by 

P\  (?) 


*i<£) 


(4.22) 


(4.23) 


(4.24) 


where 


<?<£)--  (k-lkA )  (ki-ika)  (4.25) 

and  fi()D  is  the  Fourier  transform  of  PiO?)  as  given  by  equation  (3.4),  which  in  this  case  is 


*(»>*•?+>»)  . 

Pl  0  <r0(zs)exp(xCrs-z)}  ’ 

accordingly,  we  have,  assuming  (4.11)  for  o-i, 

PiQc)  “  -  /  \  **  ^ (_^zo)  expi-ili'x'o) 


(4.26) 


(4.27) 


and 


.  ,  Eoo^^zj)  ,  r  x  ,  ,  ,  exp(~ik-xl) 

<t>\(k)  ■  I  o  -  r  r  W'fco)  «7» (-x*o) 


(2w)V0Us) 


(kt~ikA){k-ikB) 


(4.28) 


Reference  to  the  boundary  conditions  (3.18)  shows  that  A  (kx,ky)  and  B(kx,ky)  are  related  to 
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4>i  through  the  matrix  expression 


A  (kx  i  ky ) 
B(k,yky) 


N 


R  (kx  ,ky  ,z$ ) 
R  (kx,kyja) 


(4.29) 


where 


R 


1 

D 


kBexp(-kBzB)  -kBexp(-kBzs) 
-kAexp(-kAzB)  kAexp(-kAzs) 


with 

D—2kl  exp(-kzs )  exp  Cotd)  sink  (fid)  , 
d  -  z$  -  zB 
and 

R0cxtkyy r)  -  J  Ikt  Q\p(ic)  expitkjZ )  dk,  , 


(4.30) 


(4.31) 

(4.32) 


(4.33) 


ir  we  make  the  separation 

LyjOc)  ■»  £p ( Tc )  +  i  ,j  Jc)  -f  CB  (Jc)  ,  (4.34) 

where  the  terms  on  the  right  side  are,  respectively,  the  terms  on  the  right  side  of  equation 
(4.18)  substituted  individually  into  expression  (4.17),  we  see  that 

(gxf)  (4.35) 

•tr  Kj  "T"  1^0 

x  lexp{-/(fci+//ca)zs}  -  exp {-t(k,+ikB)zB)]  A  (kx,ky) 


and 


&(£)-/ 


fro(?s)eXp(r-kzs) 
2 ir 


(£*i) 


kj  -  (kg 

k,  +  ikA 


X  lexp{-i(^+^)z5)  -  cxp(-l(kt+lkA)zB}]  B (kx ,ky ) 


(4.36) 


The  remaining  integrals  are  easily  evaluated  by  contour  integration,  giving 


Cp(k)  -  -^y-  k‘V  (£x  ?)  expl-i^xa+A^yo)}  W(zo ) 


\exp(-iksz0 )  _ 

(A,  +  ikA)(k,  +  ikB) 


(4.37) 


Jl 
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(kt  -  ikA)exy(-i(k,  +  ikB)zs) 
#(*,  +  ikB) 

(kf  -  ikg)cxp(-i(k,  +  ikA)zB ) 
3/9  (k,  +  ikA) 


exp  ( -kBz0 ) 
expi-k^zo) 


CA  Qc)  -  *-lc  Ocxi)  2^^  l«xp(-/  (k,+ikB)zs)  -  exp{-/(Ai  +  ^>2^4338) 


and 


x  expl-Z^xo+Ar.^o)}  W'Cro)  Y4(kxXxo)  , 


2ir  f**  3/3 


x  (exp{-/ (fcj+M:^  )*$}  -  exp{-/(^  +  1^)2*}] 


(4.39) 


in 


x  exp(-/(AfJ(x0rj-^>>o)}  H'fco)  YB(kXfkyxo)  * 


where 


Kf(A*  ♦*,*>) 


Ar^  exp  (~k4  zs  )exp  i-kBz0) 

■  /cfl  exp  (-fcpr*  )«p  (-/f^  2q) 


(4.40) 


We  now  integrate  the  terms  of  £r  derived  above  against  o^OcJ)  over  all  within  the 
ocean  to  get  an  expression  for  Cy  for  arbitrary  erj.  If  we  define  the  transforms 


7>(F)  -  — J  dx  j  J  *  Wig)  exp (—iic-!R)  , 


(2w)3  L 


Ta  (kx,ky)  -  J  dx  J <fy  J  dz  l  <r,C3?)  K^) 


(4.41) 


(4.42) 


I 


1 


and 


x  exp{-/(Ar*x+fc,y))  exp(-k4z  )  , 


10  op  M 

r«<fc4)  -  ^  ^  *  1  <r,C?)  ^ 


(4.43) 


£aaflaafiaaaS!aaaaafl^^ 


86 


x  cxji{-l(kxx+kyy))  exp(-kBz)]  , 
we  find  that  the  equivalents  of  (4.14),  (4.37),  (4.38),  and  (4.39)  are 
Ln<£)-lOix*)-ApTr<&  , 

C0?)-£()i*0?xJ)l  vr'®- 


(fr,  +  /A^)(A*  + 
(Ai  -  ikA)oxp[-l(k,  4-  tkB)zs) 


mkt  +  ikB) 


TB(kxA) 


,  (k,- lkB)txf[-t(kt  + lkA)zB)  „  ,.s 

+  V(k,  +  ikA )  A 

(Tx*)  + 


.  x  [expH  (*,+/**  )zs}  -  exp{-/(/c,  +  lkB)zB)]  ZA  (kx,ky)  , 


(4.44) 

(4.45) 


(4.46) 


and 

fi®-**Tffx  D^rgyl  «•«> 

x  [exp{-/(^+/A^)zs}  -  exp{--/ (Ar,  4-  ^)za)l  ZB(kxtky)  , 


where 


(A'r  ,^i ) 

(A*  iky ) 


kAexp(r-kAZs )  TB(kxtky) 
kBexp(—kBzB)  Ta  (kx,ky) 


(4.48) 


Note  that  7>  Is  just  an  inverse  Fourier  transform. 


5.  Two  different  methods  for  *  general  solution 


5.1  The  solution  in  rectilinear  coordinates 

This  method  is  based  on  the  separation  of  variables  in  a  rectilinear  coordinate  system, 
and  its  final  stages  make  use  of  partial  Fourier  transforms.  The  final  result  is  presented  in  the 
form  of  a  partial  Fourier  transform,  and  is  somewhat  easier  to  convert  into  a  full  Fourier 
transform  than  into  the  untransformed  equivalent;  for  this  reason,  this  method  is  recom¬ 
mended  over  the  other  if  one  wishes  to  calculate  the  magnetic  effects  associated  with  a  front,  as 
it  interfaces  more  cleanly  into  the  mathematics  of  section  2.2  than  does  the  other  method. 

We  start  with  the  homogeneous  equation  for  any  order, 

V<0CT)  +  /<z)  ,  (5.1) 

and  assume  that  separates  as 

*(*)-  K(z) //(*>)  .  (5.2) 

Some  work  then  gives  the  equations 

kw-o  (in 

and 


VhH(xy)  +  k1  Htoy)  -  0  ,  (5.4) 

where  k 1  and  F(z)  are  an  eigenvalue  and  the  associated  eigenfunction  of  equation  (5.3)  under 
the  boundary  conditions 


dViz)  dV(z)  - 

dz  dz  r*s 


(5.5) 


derived  from  the  criteria  that  for  all  orders  2  V<f>  must  be  zero  at  the  ocean  bottom  at  zB  and 
the  ocean  surface  at  zs .  For  convenience  of  notation  we  label  the  nth  eigenvalue  and  its  associ¬ 
ated  functions  Vn(z ),  and  H„  U  Vv).  If  orthogonality  can  be  demonstrated  for  the  Vn ,  then 
they  may  be  used  in  the  construction  of  the  final  solution;  that  the  other  requirement  of  com¬ 
pleteness  is  satisfied  is  suggested  by  the  observation  that  equation  (5.1)  under  the  boundary 
conditions  (5.5)  -should  have  an  infinite  number  of  discrete  soiution3,  which  should  form  a 
complete  set  if  they  are  orthogonal  to  each  other. 
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To  demonstrate  the  mutual  orthogonality  of  the  V„  (z),  we  first  note  that,  according  to 
the  definition  of  f(z)  in  section  2.1, 


d(r0  Or) 


“/(z)  croCr)  • 


Letting  a  prime  denote  differentiation  with  respect  to  z,  we  observe  that 


~  O'mcroK,')-  VmW  +  fVmcr 0K„'  +  KMor0K„" 

~V<r0K/  +  Kmcro(n"  +  /K,') 

•  *  -  VmtroK 

by  (S.3);  given  this,  we  then  have  that 

“  <  WV-  K„VroK.)  -  (At,2  -  *m2)  Kmcr0^  , 
and  hence 


(tf-*,2)}  *V(z)v„(z)]  IJ-0  (5.9) 


by  (5.5).  Therefore,  we  have  that 

}<r0<t)Vm(*)V„Q')dz-8im  1„  , 

where  8-.  is  in  this  case  the  Kronecker  delta,  and 


h  ■  ,f  <r0(z)K)2(r)  *  . 


(5.10) 


(5.11) 


We  now  turn  to  the  nonhomogeneous  equation, 
V^OO  +  Zfe)  ^|^--p0?)  . 


(5.12) 


If  we  define  the  partial  Fourier  transforms 


I 


w  m 

0 0?)-  /  <tfr,  J  dky  F(kx>kyj)  txp{l(kxx+kyy)\ 


(5.13) 


(5.14) 


P 0?)  “  /  dkx  f  dky  P(kx,kyj)  txpil(kxx+kyy))  , 

— ~oa 

then  equation  (5.12)  is  transformed  into 

0+/<*>  -**/?-/>  ,  (5.15) 

where 

k*  -  kf  +  k*  ;  (5,16) 

equation  (5.15)  may  he  solved  through  the  construction  and  use  of  a  Green’s  function  G  such 
that 

FOctJcyj)- J  G<*j\k)  P(kx>ky*')'k'  .  (5.17) 

The  appropriate  equation  for  G  is 

^  fl(»/^)  +  /(«)  A  ^  G (:/,«  -  !H')  ;  (5.18) 

following  a  line  of  reasoning  similar  to  that  used  to  demonstrate  the  mutual  orthogonality  of 
the  VH ,  we  observe  that 

■fa  (Y*<ro G'-  Y,/<roG}  •»  V„(z)  cr0(z)  <?(z,z\Ac)  (A:2  -  A :})  (5.19) 

-  V„(f)  cro(r)8(z-z')  , 

where  a  prime  again  denotes  differentiation  with  respect  to  z,  and  hence 

0- (Ac2- Af„2)  jf  <r0(z)l^(r)G(z^',Ac)  dz  +  o>0(r')Krt  (z')  ,  (5.20) 

which  has  the  solution 

G(zj\k)  -  ■  [  kl)  cr^V'iz')  +  C(zj\k)  ,  (5.21) 

where  C  is  any  function  orthogonal  to  (?).  On  the  basis  of  symmetry  we  guess  that 
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g (jj\k)  - 1  (kj !_  k2)  \  }  y»  <*'>  . 

and  test  this  by  substituting  it  back  into  equation  (5.18)  to  get 

*-0  y« 


(5.22) 


(5.23) 


this  is  in  fact  the  expression  that  one  gets  by  using  (5.10)  to  do  a  series  expansion  of  8  Or-*') 
in  the  Vn  (z),  indicating  that  (5.22)  is  the  correct  expression  for  G. 

Accordingly,  we  have  from  (5.17)  that 


at  y  (r) 

1  (kf-k2)l  &<**,*,) 


(5.24) 


where 


i  (kx  ,ky  )  ■■  J 


<r0(z)  V„  (z )P (A*  ,kyfi)  dz  , 


(5.25) 


If  4>(£)  is  the  full  Fourier  transform  of  ^(7),  then  we  have  that 

•• 

4>(i P)-/  F(kx,kyiz)  dr  , 


(5.26) 


the  evaluation  of  which  amounts  to  replacing  the  (r)  in  (5.24)  with  their  Fourier  transforms. 
If  we  calculate^  Of)  from  (5.24)  by  a  two-dimensional  inverse  Fourier  transform  of  F,  we  get 


0Cr)-i  Mi  *„(<,,)  . 

(1-0 


MM 

A*(x^)- j*  «*c’  /  4v'J  dr'o-0(r')Kw(rU«(x'-xiy->')pOr> 

X| 

f  dk  f  dk  w\KkxX+k,y)\ 

L*'*'y)  J  *  -I  '  A„2  -  ^  -  V 


(5.27) 


,  (5.28) 


-  M(k„x,k„y)  , 


(5.29) 


where 


Mhj.)  -  t  <*>  f  *>  , 

(2tr)2  l-w2-v2 


(5.30) 
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which  reduces,  with  transformation  to  polar  coordinates  followed  by  integration  over  all  angles, 
to 

1  m 

M U&)  -  |  Mkf>)  dk  ,  (5.31) 

where 

pt-x^  +  y*  .  (5.32) 

Making  use  of  the  identities 

-  -  (-1)"  ff„(5)(z)  (5.33) 

#,<»(-*)--  (-1)"  Hnw(x)  ,  (5.34) 

where 

HHw<x)-Jn(z)+lNH(x)  (5.35) 

#«>(*),-  N*(t)  ,  (5.36) 

we  find  that 

jf  ^±^Hp(kp)dkm.2^j^^J,(kp)dk  t  (5i37) 

giving  from  equation  (5.31) 

A#0r^)  -  HP  Okp)  dk  ;  (5.38) 

this  muy  be  evaluated  by  contour  integration,  closing  the  contour  in  the  lower  half  plane  and 
excluding  the  pole  at  the  origin  as  a  nonphysical  mathematical  artifact,  to  give 

Mlxy)  -  -1/4  No(p)  .  (5.39) 

The  calculation  of  <£C?)  from  F  is  somewhat  more  tedious  than  the  calculation  of  d >(£), 
involving  two-dimensional  convolution  integrals  as  opposed  to  one-dimensional  Fourier 
transforms. 


92 


5.2  The  solution  in  cylindrical  coordinates 

This  method  is  based  on  separation  of  variables  in  a  cylindrical  coordinate  system,  and 
makes  heavy  use  of  Bessel  transforms.  For  the  final  result,  <£0?)  is  presented,  in  cylindrical 
Coordinates,  as  a  series  of  Bessel  functions  of  the  second  kind.  Disadvantages  of  this  method 
are  that  it  is  practically  suited  only  for  electrical  conductivity  perturbation  distributions  localized 
about  the  origin,  since  an  extended  distribution  requires  more  terms  in  the  series  expansion 
than  does  a  localized  one,  that  large  numbers  of  terms  are  necessary  to  adequately  describe  the 
electric  field  close  to  the  origin,  since  successively  higher  order  terms  in  the  expansion  for  <f> 
tend  to  be  progressively  more  concerned  with  the  field  in  the  neighborhood  of  the  origin  and 
less  with  the  field  at  great  distances,  and  that  using  this  method  to  calculate  magnetic  effects  is 
more  involved  than  using  the  method  discussed  in  the  last  section;  balanced  against  this  is  the 
availability  of  numerically  swift  algorithms  for  computing  Bessel  transforms  (Anderson,  1979; 
Chave,  1983).  Accordingly,  this  method  is  at  its  best  in  computing  the  electric  field  associated 
with  a  localized  electrical  conductivity  variation,  particularly  the  field  at  large  distances  from  this 
variation. 

We  start  by  assuming  the  separation 

.  ,  (5.40) 

which,  when  substituted  into  the  homogeneous  equation  (5.1)  yields  for  V  and  H 

+  /(z)  ~~  -  k2  f'(z)  -  0  (5.41) 

and 

<p,0)  +  k1  Hipfi)  -  0  ;  (5.42) 

the  boundary  conditions  on  <£,  that  2'V<A-0  at  the  ocean  bottom  at  zg  and  the  ocean  surface  at 
zs ,  give  that 

dV(z) 
dz 


iki 
<k  n 


(5.43) 


serving  to  restrict  solutions  of  (5.41)  to  discrete  eigenvalues  k%  and  associated  eigenvalues 
K„(z).  We  observe  that  (5.41)  and  (5.43)  are  identical  in  form  to  (5.3)  and  (5.5)  respectively 
in  the  previous  section,  Indicating  that  the  eigenfunction  orthogonality  arguments  made  there 
can  be  extended  to  the  work  of  this  section. 


We  then  define  the  transforms 

eo  ** 

£  f  kdk  Fm(z,k)elm»Jm(kp) 

m- 0  0 


end 


(5.44) 


pip  fit)  -  £  J  k  dk  Pm(z,k)  elmt  Jm(kp ) 

m- 0  0 


(5.45) 


«nd  their  inverse  transforms 

Fm(z,k)--±}apdp^dO<t>(pfij)e-‘"'jm(lcp)  (5.46) 

and 

PmiXsk)  -  |  p  dp  |  d*  pip fij)  4  (*p)  ,  (5.47) 

where  ^  andp  are  as  in  equation  (5.12).  Substituting  expressions  (5.44)  and  (5.45)  into  equa¬ 
tion  (5.12)  and  using  the  identity 

Vhel”'tJm(kp)--k2e,m>Jm{kp)  (5.48) 

gives  the  equation 

-£iFm(x,k)  +  f(,)-jLFm<z,k)-kiFmiz,k)-Pm(z,k)  ,  (5.49) 


which  is  seen  to  have  the  same  form  us  equation  (5.15)  in  the  previous  section. 

Because  of  the  observed  similarities  between  the  work  of  this  section  and  that  of  seo 
tion  4.1,  we  can  immediately  write 


G(zj',k)  Pm(z\k  dz') 


(5.50) 


where  G  is  as  given  by  equation  (5.22),  and  hence 

I  -Tir-fT  —j—  K*,Ak) 
»- 0  l» 


(5.51) 


where 


Substituting  (5.51)  back  into  the  transform  (5.44)  then  gives 


oo  oo  V  (z) 

*^)-H  -y-1  e‘m 6  (5.53) 

mm0  rt-mO  '» 

where 

I^lp)-/ kdk^f~Jm{kp)  .  (5.54) 

o  **  —  * 

We  can  use  the  identities  (5.33)  through  (5.36)  in  section  4.1  together  with  the  observation 
that 


Km(-k)  -  (-l)m  Km(k)  (5.55) 

to  demonstrate  that 

Lm  (p)  -  #  2  k  dk  jfrj?  H”)(kp)  ’  (5-56) 

which  may  be  evaluated  by  contour  integration,  closing  the  contour  in  the  lower  half  plane  and 
detouring  the  section  along  the  real  axis  below  the  origin,  to  get 

L<m  (p)  •*  —  ir  Km  (k„)  Nm  (k„p)  .  (5.57) 

Note  that  each  of  the  Km  need  be  evaluated  for  only  a  single  value  of  k  in  these  calculations. 
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